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I.  INTRODUCTION 


\  - 

V 

The  normal  mode  model,  NEMESIS,  was  designed  and  implemented  at  the 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin  (ARL:UT), 
to  aid  in  the  investigation  of  low  frequency,  range  invariant,  acoustic 
propagation  effects.  The  mathematical  basis  and  numerical  scheme  of  this 
model  are  described  by  Gonzalez  and  Hawker  (1980) .  NEMESIS  computes  the 
eigenvalues  and  normal  modes  of  the  depth  separated  acoustic  wave  equation. 
Group  velocities  and  modal  attenuation  terms  are  also  computed.  Propaga¬ 
tion  loss  and  velocity  potential  for  a  source  and  receiver  at  a  given 
range  can  be  calculated  from  these  quantities  (Gonzalez  and  Hawker,  1980). 
A  computer  model,  PLMODE,  was  designed  to  perform  these  computations,  and 
the  computer  code  for  both  NEMESIS  and  PIMODE  was  prepared  for  export  to 
other  research  laboratories.  This  report  is  designed  to  accompany  this 
software  and  to  assist  persons  who  wish  to  use  it.  The  contents  include 
descriptions  of  the  modeled  environment,  the  input  data  and  their  format, 
and  the  output  produced  by  the  models.  An  explanation  of  the  diagnostic 
messages  that  may  arise  during  execution  of  the  models  is  also  included. 
Estimates  of  memory  requirements  and  execution  times  are  given  along 
with  sample  data  sets  and  their  resulting  output. 

r\ 

Although  it  provides  information  beyond  the  scope  of  a  user's  manual, 
an  appendix,  "Installation  Information,"  has  been  added  to  this  document 
to  aid  in  the  retrieval  of  the  software  from  the  magnetic  tape  on  which 
it  was  written  for  export. 
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.  II.  NEMESIS:  DOCUMENTATION 


The  following  chapter  is  extracted  from  the  computer  model  NEMESIS. 


I 


II-l 


NEMESIS  ARL'.UT  normal  MOOE  model 
VERSION  Tt2  FEBRUARY  i960 

DESIGNED  AND  IMPLEMENtED  BY 

RUTH  GONZALEZ 

DOCUMENTATION  for  nemesis  I*  GIVEN  IN 
THE  ACOUSTIC  NORMAL  MODE  MODFL  NEMESIS 

BY 

RUTH  GONZALEZ  AND  KENNETH  F.  HAWKER 
ARL'.UT  TECHNICAL  REPORT  ARL-TR-80-13 
AND 

user's  MANUAL  FOR  NEMESIS  And  PLMODE 
BY 

RUTH  GONZALEZ  AND  SUSAN  G.  PAYNE 

arl'.ut  technical  memorandum  apl-tm-bo-6 
the  APPLIED  RESEARCH  LABORATORIES 
the  university  of  Texas  at  Austin 

SPONSOREO  BY  NOROA  CODF  5 ?«*  NSTL  STATION*  MISSISSIPPI  39529 


MATHEMATICAL  and  physical  DESCRIPTION 


nemesis  is  low  frequency*  range-invariant*  wave  theory  normal  mode 
model*  it  WAS  DESIGNED  TO  COMPUTE  EIGENVAl  UeS  ano  normal  mooes  in  a 
horizontally  stratified  DEEP  ocean  with  single  channel  profiles  ano 

MULTIPLE  FLUID  SEDIMENT  LAYERS  OVERLYING  A  SUBSTRATE.  SOUND  SPEEDS 
IN  THE  WATER  AND  SEDIMENT  VARY  WITH  DEPTH.  THE  DENSITY  IS  CONSTANT 
UITHIN  EACH  LAYER*  THE  LAST  LAYER  IS  A  HOMOGENEOUS*  SEMI-INFINITE* 
FLUIO  OR  SOLID  SUBSTRATE*  I*E.»  THE  COMPRERSIONAL  ANO  SHEAR  SPEEDS 
AND  DENSITY  REMAIN  CONSTANT  WITH  DEPTH.  SnUND  SPEEDS  AND  DENSITIES 
CAN  BE  DISCONTINUOUS  at  layer  INTERFACES,  the  model  NUMERICALLY 
Computes  the  oiscrete  eigenvalues  and  eigenfunctions*  or  normal  modes* 
OF  THE  DEPTH  SEPARaTEO  ACOUSTIC  WAVE  EQUATION  FOR  THIS  TYPE  OF 
ENVIRONMENT*  THE  CONTINUOUS  PART  OF  THE  SPECTRUM  OF  THE  SOLUTION 
IS  ignored*  group  velocities  are  ALSO  COMPUTED* 


DEPTH  EQUATION 


2  2  2  2 
(D  /  DZ  1  U  ♦  (K  <Z)  -  K  )  U  *  0 
N  N  N 


WHERE  U  IS  THE  EIGENFUNCTION 
N 

K  IS  THE  EIGENVALUE 
N 

K (Z)  ■  2*PI*F  /  C<Z) 

c<z>  is  the  sound  speed 
F  IS  THE  source  FREQUENCY 


ANY  QUESTIONS  CONCERNING  ThE  APPLICABILITY  OF  THE  MODEL  TO  A 

particular  situation  should  be  directed  to 


KENNETH  E.  HAWKER 

APPLIED  RESEARCH  LABORATORIES 

THE  UNIVERSITY  OF  TEXAS  AT  AUSTIN 

P*  O*  BOX  8029 

AUSTIN*  TEXAS  T87]2 

TEL*  (512)  835-3AQ7  OR  835-3200 


ON  INPUT 


The  FOLLOWING  RESTRICTIONS  are  IMPOSED  On  the  INPUT  to  nemesis 

(1)  there  must  be  at  least  ?  layers  (water  and  substrate)  but  no 
more  than  io. 

(2)  ALL  DEPTHS  MUST  BE  IN  INCREASING  OROER. 

O)  depths  at  later  inteRPaccs  must  be  pqual, i*e.»  the  last  depth 
at  THE  bottom  of  a  LAYER  MUST  EQUAL  ThE  FIRST  DEPTH  AT  TOP  OF 
THE  NEXT  LAYER. 

(4)  FIRST  SPECIFIED  DEPTH  MUST  BE  0  (AIR-SURFACE  INTERFACE). 

(5)  THE  SOUNO  SPEED  MINIMUM  IN  The  WATF»  and  sediments  MUST  BE  less 
than  the  SHEAR  ANO  cOMPRtSSlONAL  WAVE  SPEEDS  IN  THE  SUBSTRaTE. 

(6>  SHEAR  SPEED  OF  ZERO  IN  THE  SUBSTRaTf  IMPLIES  A  FLUID  SUBSTRATE 
AND  A  NONZERO  SHEAR  SPEED  IMPLIES  A  SOLID  SUBSTRATE.  NO  SHEAR 
SPEEDS  ARE!  ALLOWED  IN  THt  WATER  Or  SEDIMENTS. 

(7)  the  shear  speed  must  be  less  than  ime  compressional  wave  spEto 

IN  THE  SUBSTRATE. 

(8)  IP  A  PRESSURE  RELEASE  BOTTOM  IS  SPECIFIED#  THERE  ARE  AN 
INFINITE  number  of  MODES.  The  COMPHFSSIONAL  speed  is 

used  AS  THE  PHASE  VELOCITY  CUTOFF,  IN  THIS  CASE.  AND  IS  NOT 
USED  IN  THE  BOUNDARY  CONDITIONS. 

(9)  FOR  SOME  OATA  INPUTS#  DEFAULTS  ARE  « IVEN  IN  PARENTHESES. 

these  can  be  invoked  by  leaving  the  appropriate  data  field 
blank  or  setting  it  to  «. 

(10>  MAXIMUM  NUMBER  OF  DEPTH-SOUND  SPEED  PAIRS  IS  125. 

(ID  MAXIMUM  MODE  wHlCH  CAN  BE  COMPUTED  TS  MODE  NUMBER  300. 

(12)  MAXIMUM  NUMBER  OF  MESH  POINTS  PER  mode  IS  2000.  THE  NUMBER  OF 
MESH  POINTS  IN  A  LAYER  IS  EQUAL  To  THf  NUMBER  OF  M£SH  INTERVALS 
IN  THAT  LAYER  PLUS  1. 

(13)  THE  NUMBER  OF  mesh  INTERVALS  in  Each  LAYER  MUST  BE  EVEN  AND 

greater  than  or  equal  to  4.  there  are  no  mesh  intervals  in  )he 

SUBSTRATE.  the  SOLUTION  is  ANALYTICALLY  COMPUTED  IN  That 

region. 

(14)  ALL  GUANTITiES  MUST  BE  SPECIFIED  in  The  MKS  SYSTEM  EXCEPT 
DENSITY,  WHICH  IS  IN  ThF  CGS  SYSTEM. 

(15)  the  DEPTHS  ASSOCIATED  WITH  BOTH  ThE  COMPRESSIONaL  AND  SHEAR 
SPEEDS  in  The  SUBSTRATE  are  the  Substrate  INTERFACE  depths. 

The  COMPRESSIONAL  SPEEo  IS  FIRST  ANn  THE  SHEAR  SPEED  IS  SECOND. 


INPUT  TO  NEMESIS  IS  VIA  FORMATTED  READ  STATEMENTS.  EACH  ITEM  MUST 
BE  POSITIONED  IN  THE  CORRECT  FIELD  AS  GIVEN  bY  THE  FORMAT  STATEMENT. 

integer  values  and  exponents  must  be  right.justified  in  their  field. 


CARD  COLUMN  NO 

NO  OF  CARO  DESCRIPTION  OF  DATA  FORMAT 


l  i-80  header  information  used  to  describe  data 


baio 


2 


3 


4 


5 


1-10 

SOURCE  FREQUENCY  IN  HZ 

Flu. 3 

11-20 

depth  where  the  integrations  are  to  be 
must  BE  in  WATER  (oefault  is  at  sound 
MINIMUM  IN  the  WATER  OR  MIDDLE  OF  The 
COLUMN.  DEFAULT  is  REcOMMFNDED) 

matched* 

SPEED 

WATER 

Flu.3 

21-30 

wave  NUMBER  k 1 

Flo. 3 

31-40 

wave  number  K2 

Flu. 3 

41-50 

wave  NUMBER  INCREMENT.  Kl*  k2»  and  INCREMENT 

Flo.3 

are  for  additional  diagnostic  information,  the 
differential  equation  is  integrated  using  the 
have  NUMBERS  K1  to  K 2  IN  The  INCREMENT  SPECIFIED 
IN  COLS  41-50.  THE  DISCONTINUITY  IN  THE  DEPTH 

derivatives  of  the  trial  modes  is  tabulated  for 
each  wave  number,  not  recommended  unless  tracing 
CONVERGENCE  OF  numerical  PROCESS.  all  three 
QUANTITIES  must  BE  POSITivf  to  activate 

diagnostics 

1-5  LOWEST  MODE  NUMBER  TO  be  COMPUTED  <D£FAULT  IS  1)  15 

6-10  HIGHEST  MODE  number  to  be  COMPUTED  (DEFAULT  IS  15 

TO  COMPUTE  all  MODES  UP  TO  MODE  NUMBER  300) 

OPTIONS  TO  BE  USED  <1*YES*  QsNO) 

5  CORRECT  SOUND  SPEEDS  FoR  CmrvATUrE  oF  EARTH  15 

10  USE  PRESSURE  RELEASE  SUBSTRATE  15 

)5  SUPPRESS  CALCULATION  OF  HlnDEN  DEPTHS*  NOT  15 

RECOMMENDED 

20  TURN  ON  DIAGNOSTIC  PRINTS  AS  MODES  ARE  COMPUTED.  15 

not  Recommended  unless  tracing  convergence  of 

NUMERICAL  PROCESS 

l-IO  DENSITY  of  LAYER  Flu. 3 

11-15  number  of  MESH  INTERVALS  In  LAYER  (DEFAULT  IS  15 

MAXIMUM  ALLOWABLE) 


II-6 
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6A  t  B*  • 


1-10 

SOUND  SPEED  PROFILE  IN  ThIS 
DEPTH 

layer 

FlO. 3 

11-20 

SOUND  SPEED 

F 1 0  •  3 

21-30 

DEPTH 

Flu. 3 

•  •  • 

•  •  • 

... 

71-80 

SOUND  SPEED 

F  i  o  .3 

1-10 

BLANK  CARD 

REPEAT  CARDS  5-7  FOR  EACH  CEDANT  LAYER 

DENSITY  OF  SUBSTRATE 

FlO. 3 

1-10 

depth  of  sediment/substRatp 

INTERFACE 

F10.3 

11-20 

COMPRESSIONAL  SPEED 

Flo  .3 

21-30 

OEPTH  OF  SEDIMENT/SUBSTRATF 

INTERFACE 

Flo  *3 

31-40 

SHEAR  SPEED 

Flo. 3 

II-7 


V 

v* 


ON  OUTPUT 


NEMESIS  USES  FILES  LFOUT,  lTAPEI*  and  lfdug  as  output  files 
ON  LFOUT 


NEMESIS  USES  FILE  LFOUT  to  OUTPUT  ALL  PRINTED  INFORMATION  such  AS 
inputs  and  ouantities  that  were  calculated  in  the  program  such  as 
EIGENVALUES*  group  velocities*  and  ATTENUATION  terms,  the 
ATTENUATION  terms  for  each  mode  are  LISTED  in  the  following  order* 

attenuation  in  water  LAYER 
attenuation  in  first  sediment  layer 

ATTENUATION  in  last  sEDIMLNT  LAYER 
ATTENUATION  in  substrate  t compress I On Al ) 

ATTENUATION  in  substrate  <shear> 

ON  LFDIAG 


nemesis  uses  LFDIAG  For  PRINTED  diagnostic  MESSAGES.  CAUTION  -  IF 
DIAGNOSTIC  PRINT  option  IS  TURNED  On*  this  FILE  CAN  BE  VERY  LARGE. 

The  FOLLOWING  ERROR  ANO  OIaGNOSTIC  MESSAGES  a RE  GENERATED  BY  NEMESIS: 

FATAL  INPUT  ERRORS 

ONLY  lo  LAYERS  ARE  ALLOWED 

at  LEAST  TWO  LAYERS  MUST  RE  USED  TO  RUN  NEMESIS 
SOUND  SPEED  ARRAY  OVERFLOW 
all  depths  MUST  be  in  INCREASING  order 
depths  at  layer  INTERFACES  MUST  be  EQuai 

each  LAYER  MUST  HAVE  AT  LEAST  TWO  DEPTh-SoUnO  SPeFD  POINTS 
depth  at  air-wateh  INTERFACE  must  be  zero 


SOUND  SPEED  MINIMUM  of  the  water  AND  SEOImEnT  LAYERS 
must  be  less  than  THE  COMPRESSIONAL  VELOCITY 

shear  speed  must  be  less  than  the  comprfssional  speed 

IF  THE  SHEAR  SPEED  IS  NONZFRO*  THE  SOljNn  SPEED  MINIMUM 
must  be  LESS  than  THE  ShEAP  SPEED 

OPTIONS  MUST  BE  SPECIFIED  wllH  A  0  OR  1 

TOO  MANY  SPEEDS  FOR  SuRSTwaTl 

FATAL  ERRORS  (NUMERICAL  AND  PROGRAMMING ) 

integration  count  exceeds  maximum  allowable*  MODE  NO  = 

ALL  MATCH  POINTS  EXHAUSTED*  CONVERGENCE  NOT  ACHIEVED*  FOR  MODE  NO 

ERROR  IN  NUMROV*  CALLED  FROM  INTGRl *  IEP  s 

ERROR  IN  SUBROUTINE  INTrP  AFTER  CALL  TO  IcSlCU 

ERROR  IN  SUBROUTINE  INTRP  AFTER  CALL  TO  IcSEVU 

ERROR  in  CALL  To  Secant  FOR  MODE  NUMBER* 

ERROR  IN  COMPUTATION  OF  RAYLt-IGH  VELOCITY 

ERRORS  DETECTED  IN  IMSL  RoUTinlS 

***  TERMINAL  error*  ier  *  from  IMSl  routine 

•••  WARNING  WITH  Fix  ERROR  ier  =  From  jmsL  ROUTINE 

***  WARNING  ERROR  IER  =  FROM  IMSL  ROUHNE 

***  undefined  error  iep  =  from  imSl  routine. 


i* 


OIiGNOSHC  PRINTS  FOR  TRACING  CONVERGENCE  OF  NUMERICAL  ACCESS 
(FOR  MORE  DETAILED  INFORMATION  SEE  GONZALfZ  AND  HAWKER*  I960) 

the  current  mode  number  being  computed* 
mooeno  * 

ONE  OF  THE  FOLLO^nG  IS  PRINTED  FOR  EACm  TRIAL  EIGENVALUE* 

FIRST  ESTIMATE 

second  estimate 
secant  iterate 

right  BOUND  ESTIMATE  from  refine 
left  bound  ESTIMATE  FROM  refine 
right  BOUND  ESTIMATE  from  zeroin 
left  bound  ESTIMATE  FROM  zeroin 

BISECTION  ESTIMATE 

THE  TRIAL  EIGENVALUE*  NUMBER  OF  ZEROS  Im  THE  TRIAL  EIGENFUNCTION* 
and  INTEGRATION  NUMBER* 


NZEKOS 


INTEGRATION  no 


THE  FIVE  POINTS  CEnTEREO  AT  THE  MATCH  PoInT  FOR  EACH  UPPER  AND 
LOWER  HALF  OF  TRIAL  MODE* 


U-UPPER 


U-LOWER 


the  depth  derivative  Evaluated  at  the  match  point  f0r  each 
UPPER  AND  LOWER  HALF  OF  TRIAL  MODE* 


UPPER  OERV 


LOWER  OERV  s 


the  DISCONTINUITY  In  THF  DEPTH  DERIVATIVES  AT  THE  match  POINT* 
I.E..  Disc  «  UPPER  UERV  -  LOWER  DERV.  RATlO  IS  THE  SCALING 
FACTOR  THAT  WAS  APPLIED  to  THE  UPPER  HAI  F  To  FORCE  CONTINUITY 
IN  The  UPPER  and  LOWER  HALVES  OF  TRIAL  mOdE* 


DISCONTINUITY 


RATIO  = 


t  i 

t 

/  i 


If  i 


11-10 


DIAGNOSTIC  PRINTS  FOR  TRACKING  DISCONTINUITY  <phi>  in  first 
DERIVATIVE  of  TRIAL  EIGENFUNCTIONS  for  a  range  of  trial 
eigenvalues  for  a  FIXED  match  points 


PHI  ZEROS  ZMATCH 


ON  LTAPE1 


nemesis  generates  a  disk/tape  file  on  ltarfi  that  contains  all  of  the 
RELEVANT  INPUT  data  TOGETHER  with  ALL  of  The  EIGENVALUES*  GROUP  VELOC¬ 
ITIES*  ATTENUATION  TERMS*  and  NORMAL  MODES.  NOTE  -  WITH  THE 
EXCEPTION  of  THE  HEADER  INFORMATION*  all  UMAnTiTIES  are  OUTPUT  AS 
FLOATING  POINT  NUMBERS.  LTaPE]  HAS  THE  FOL» OwlNG  STRUCTURE* 


LOGICAL  WORD  no  DESCRIPTION  OF  RECORD  NO  OF  WUkDS 

RECORD  NO  IN  RECORD  PER  RECORU 


1 

1-8 

header  information 

8 

2 

1 

NUM8£R  OF  LAYERS 

20 

2 

maximum  allowable  layers 

3 

3.0 

4 

frequency  in  hertz 

5 

COMPRESSIONAL  SPEED  IN  SUBSTraTe 

6 

shear  speed  in  substrate 

7 

Rayleigh  spled  at  substrate  interface 

8 

pressure  release  substrate  <isyes.o*no> 

9 

last  modl  COMPUTED 

10 

FIRST  MODE  COMPUTED 

11-20 

o.o 

3 

1-10 

INDEX  D F  THt  LAST  MESH  POINT  IN  EACH  LAYER 

of  the  modes,  beginning  with  layer  i 

NO 

4 

1-10 

INDEX  OF  THt.  LAST  POINT  IN  EACH  LAYER  OF 

The  INPUT  SOUND  speed  PROFILE 

NO 

5 

1-10 

depth  at  each  layer  interface 

NO 

6 

1-10 

Mesh  InCPFMcNT  (In  ME1FRS)  FOR  EaCh  LAVER 

NO 

7 

1-10 

density  of  lach  layer 

NO 

8 

1-N1 

DEPTHS  CORRESPONDING  TO  THE  INPUT  SOUND 
speed  profile 

Nl 

9 

1-N1 

SOUNO  SPEEDS  CORRESPONDING  TO  T«E  INPUT 

souno  speed  profile 

Nl 

10  - 
(N3*9) 

1-N2 

NORMAL  MODES  EVALUATED  AT  EACH  MESH  POINT 
(THERE  IS  ONE  RECORD  PFR  MODE) 

N2 

N3*10 

1-N3 

eigenvalues  corresponding  to  each  mode 

N3 

N3*ll 

1-N3 

GROUP  VELOCITIES  OF  EACH  MODE 

N3 

N3M2 

1-NS 

attenuation  terms  per  i  ayEr  per  mode 

attenuation  terms  for  fach  mode  in  each 

NO 

LAYER,  BEGINNING  WITH  I  AYER  1*  BOTH  SHEAR 

and  compressional  attemuaTions  are  COMPUTED. 
therefore,  there  are  a*  many  attfnuation 
terms  for  each  mode  as  there  are  layers*i 

(N5*l ) -2N5  ATTENUATION  TERMS  FOr  «?EcOnD  MOOE  COMPUTED 
<2N5*1)-3N5  ATTENUATION  TERMS  FOR  THIRD  MODE  COMPUTED 


(  <N3-1)*N5*1)  ATTENUATION  TERMS  FOR  i  AST  MODE  COMPUTED 
-N3*N5 


end  -  of  -  file 


TOTAL  NUMBER  OF  WORDS  IN  FILL  LTAPE1  IS 
28  ♦  5*N0  ♦  2*Nl  ♦  (2*N?*N5) *N3 


DEFINITIONS 

NO  MAXIMUM  ALLOWABLE  LAYERS 

N1  TOTAL  NO.  OF  POINTS  IN  INPUT  SOUnO  SPEED  PROFILE.  ACCESS 
WORD  NUMBER  (LaYERS-1)  IN  ATh  RECOHo 


II-I2 


TOTAL  NO*  OF  MESH  POINTS  IN  MOOES.  ACCESS  WORD  NUMBER 
(LAYERS-!)  IN  3RD  RECORD 

TOTAL  NO*  OF  MOOES.  SUBTRACT  FIRST  MODE  COMPUTED  FROM  LAST 
MODE  COMPUTEO  AND  AO»  1 

N3  •  NUMBER  OF  LAYERS 

number  OF  LAYERS  ♦  1 


REQUIRED  RESIDENT  ROUTiNFS 


STANDARD  FORTRAN 

ABSt AL0610* AMAXl • AMlNlf FLOAT tIABS»lFlX,MAXO>MINOtMODtSQRT 


NONSTANDARD  -  THESE  RUST  BE  WRITTEN  AT  EACm  nEW  INSTALLATION 
EOF  (U> 

purpose  function  that  tests  for  an  end  of  file 

CONDITION  FOLLOWING  A  READ 
FORMAL  PARAMETERS  U  -  UNIT  (FILE  NAME) 

VALUE  RETURNED  ZERO  if  NO  END  OF  FILE  IS  ENCOUNTERED 

nonzero  otherwise 

IBUFOUT  (U»IFRhAT* ARRAY »NPTS) 

PURPOSE  FUNCTION  THAT  BUFFERS  OUT  A  CONTIGUOUS 

BLOCK  OF  WORDS 

formal  parameters  U  -  UNIT  (FIlE  name) 

Ifrmat  -  ONE  For  BINARY  format 

ARRAY  -  ADDRESS  oF  FIRST  WoRD  IN  BLOCK 

NPTS  -  NUMBER  OF  WORDS  TO  BE  BUFFERED  OUT 

value  returned  zero  if  SUCCESSFUL 

NONZERO  OTHERWISE 


■;xv  I 


CURRENT  LIMITS 


VARIABLE  DESCRIPTION 

meshmx  maximum  number  of  mesh  Points 

lyrmax  maximum  number  of  layers 

mxinpt  maximum  number  of  DEPTH-SOUND  speed  PAIRS 

WHICH  CAN  BE  INPUT 

NMOOMX  MAXIMUM  mode  number  which  CAN  BE  COMPUTED 


LIMIT 


200L 

ly 

1Z& 


ARRAY  DIMENSIONS 


ARRAY 

AKSQR 

DEPTH1 

speed 

U 

AKSURH 

BUFFI 

BUFF2 

BUFF3 

CHALF 

dElz 

H 

MESH 

NET 

NPTS 

RHO 

atten 


COfcF 

bufakn 

bufgv 

bufatt 

iheadr 
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Array  dimension 

MESHMX 

MESHMX 

MeShMX 

MeShMX 

LYRmAX 

LYRmAX 

LyRmAX 

LyRmAX 

LYRMAX 

LYRmAX 

LYRMAX 

LYRmAX 

LYRmAX 

LYRmAX 

LYRmAX 

I YRMAX*1 

MXINPT 

MXINPT 

MX  I NPT*i 

NmODMX 

nmOomx 

NMODMX* (lYRMAX^I) 

bo  characters 


•V-  <,  'i , 


COMPUTER  RESOURCES 


NEMESIS  MAS  programmed  On  a  CONTROL  DATA  CORPORATION  cyber  171 • 

UNDER  The  nOS  1,3  LEVEL  48S  OPERATING  SYSTEM.  THIS  COMPUTER  HAS 
A  60-BIT  WORO  WITH  10  CHARACTERS  PER  WORD.  THE  FORTRAN  COMPILER 

is  fortran  four  extended  u,7.4B5>  which  is  described  in  cdc 
publication  604978oo»  revision  d.  efforts  were  made  to  maintain 

ANSI  STANDARD  FORTRAN  CODE.  PORTIONS  OF  The  CODE  ARE  UNAVOIDABLY 

Installation  dependent*  such  as  the  number  of  characters  per  word. 

ON  THIS  MACHINE*  THE  last  WORD  ADDRESS  (LWA)  of  nemesis  was 
62121  OCTAL  WORDS,  ThIS  CAN  BE  REDUCED  BY  DECREASING  THE  SIZES 
OF  THE  ARRAYS  IN  ThE  MAIN  PROGRAM  AND  APPROPRIATE  COMMON  BLOCKS* 

as  described  in  the  section*  current  limits. 

nemesis  was  run  with  a  typical  deep  water  profile  with  two 
sediment  layers  and  a  fluid  substrate,  eigenvalues,  eigenfunctions. 
Group  velocities*  aNO  attenuation  terms  were  COMPUTED  and  stored  on  a 
FILE.  the  total  EXECUTION  time  for  the  140  MODES  with  1999  mesh 
POINTS  PER  MODE  was  194.15  CP  StCONDS  (OPT TMlZED) •  SEE  ARL-TM-BO-6* 
CHAPTER  3.  FOR  FURTHER  DETAILS  OF  THIS  EXAMPLE* 


III.  NEMESIS:  EXAMPLES 


Two  data  sets  and  the  corresponding  printed  output  from  NEMESIS  are 
included  in  this  chapter.  The  modeled  environment  is  shown  in  Fig.  1. 

It  consists  of  four  layers:  water,  two  sediment  layers,  and  the  substrate. 

3 

The  density  in  each  of  the  four  layers  is  1.0,  1.5,  1.75,  and  2.0  g/cm  , 
respectively.  The  sound  speeds  in  the  sediment  layers  vary  linearly  with 
depth,  and  there  is  a  discontinuity  in  sound  speed  at  the  water-sediment 
interface.  Note  that  the  sound  speed  must  be  specified  at  the  top  and 
bottom  interface  of-  each  water  and  sediment  layer,  and  that  within  the 
substrate  the  shear  and  compressional  speeds  must  be  constant. 

The  first  data  set  does  not  include  shear  waves  in  the  substrate; 
i.e.,  the  substrate  is  a  fluid,  and  thus  the  shear  speed  is  zero.  At 
50  Hz  there  are  149  modes,  as  shown  in  the  printed  output.  The  five 
attenuation  terms  listed  for  each  mode  pertain  to  the  compressional  wave 
in  the  water,  the  sediment  layers  and  the  substrate,  and  the  shear  wave 
in  the  substrate,  in  that  order.  Because  shear  waves  are  not  included, 
the  shear  wave  attenuation  term  is  identically  zero  for  all  modes.  In 
this  data  set  the  option  to  suppress  hidden  depth  calculations  is  not 
chosen.  Thus,  for  modes  with  turning  points  high  above  the  water-sediment 
interface  (modes  1  to  49),  the  sediment  attenuation  terms  are  zero.  Note 
that  the  phase  velocities  of  the  modes  are  bounded  by  the  minimum  sound 
speed  in  the  water  and  the  substrate  compressional  speed.  This  computa¬ 
tion  requires  194  cp  sec  on  a  Control  Data  Corporation  CYBER  171  when 
NEMESIS  is  compiled  under  maximum  optimization. 

The  second  data  set  describes  the  environment  in  Fig.  1  with  a  solid 
substrate.  The  substrate  compressional  speed  is  5000  m/sec  and  the  shear 
speed  is  2700  m/sec.  At  50  Hz  there  are  282  modes  with  phase  velocities 
bracketed  between  the  minimum  sound  speed  in  the  water  and  the  shear 
speed.  This  computation  requires  347  cp  sec  on  the  CYBER  171  in  the 
optimized  mode. 


-*-%  m  'f  i**,-'*’- '  >X**~  'fi 64NV  i  >«-n 


SHEAR  COMPRESSIONAL 
SPEED  SPEED 


FIGURE  1 

A  DIAGRAMMATIC  REPRESENTATION  OF  THE  ENVIRONMENT 
DESCRIBED  IN  THE  DATA  SETS  FOR  NEMESIS 


ARL:UT 
AS40-1066 
SOP. GA 
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institutes  the  data  set  of  the  first  example 
[el  NEMESIS,  for  this  report. 
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Ml'NK  PROF  I  LI ,  Cl  -  1500,  Z1  -  1000,  IPS 
50. 


0074,  S  -  1233.33 


1.000 

0.000  1527. 061  51.000  1523.539 

202.000  1515.024  253.000  1512.728 

404.000  1507.351  455.000  1505.952 

606.000  1502.836  657.000  1502.085 

808.000  1500.598  859.000  1500.314 

1010.000  1500.001  1061.000  1500.053 

1212.000  1500.587  1263.000  1500.880 

1414.000  1502.024  1465.000  1502.492 

1616.300  1504.076  1667. 00e  1504.669 

1818.000  1506.570  1869.000  1507.254 

2020.00C  1509.363  2071. ee0  1510.133 

2222.000  1512.426  2273.000  1513.223 

2424.000  1515.635  2475.000  1516.465 

2626. oee  1518.963  2677.000  1519.618 

2628.800  1522.377  2879.000  1523.249 

3030.000  1525.853  3081.000  1526.738 

3232. ee?  1529.374  3283. tee  1530.266 

3434.000  1532.926  3485.000  1533.827 

3636.000  1536.503  3667.000  1537.408 

3636.000  1540.095  3689.000  1541.005 

4040.00  0  1543.700  4091  .-000  1  544.612 

4242.000  1547.314  4293.000  1548.227 

4444.000  1550.934  4495.000  1551.849 

4646.000  1554.558  5000.000  1560.917 

1 .500 

5000.000  1550. ee?  5100.000  1550.000 

1.750 

5100.000  1650.000  5200.000  1680.000 


101.000  1520.411  152.000  1517.543 

303.000  1510.725  354.000  1508.915 

505. 00e  1504.760  556.000  1503.712 

707.000  1501.477  758.000  1500.978 

909. 00e  1500.127  960.000  1500.024 

1111.000  1500.170  1162.000  1500.352 

1313. 0ee  1501.216  1364.000  1501.603 

1515.000  1502.985  1566.000  1503.521 

1717.000  1505.276  1766.000  1505.919 

1919.000  1507.943  1970.000  1609.662 

2121.000  15ie.S80  21”2.000  1511.655 

2323.000  1514.013  2374.000  1514.828 

2525.000  1517.296  2576.000  1518.130 

2727.000  1520.561  277 e.ae0  1521.525 

2929.000  1524.108  2980.000  1524.988 

3131.000  1527.608  3182.000  1528.499 

3333.000  1531.147  3384. e00  1532.045 

3535.000  1534.712  3586.000  1535.616 

3737.008  1536.297  3786.000  1539.205 

3939.000  1541.897  3990.000  1542. 8e7 

4141.000  1545.506  4192.000  1546.419 

4343.000  1549.123  4394.000  1550.037 

4545.000  1552.746  4596.000  1553.661 


HI  PIS  IS  ARL-UT  NORMAL  MODI  MODEL 
TIBS  ION  1.2 


MUNi  profile,  ci  »  1500.  zi  «  1000,  ips  *  .00?*,  b  *  1233.33 


FREQUENCY  (HZ)  50.000 

NL'MBIB  OP  EIISTINC  MODIS  146 

FIRST  MODI  COMPUTED  1 

LAST  MOEE  COMPUTED  149 

DEPTH  (M)  AT  WHICH  MODES  ARE  MATCHID  1010.000 


NO  OPTIONS  SELECTED 


GiOACOUST I C  PROEILE 


SHEAR  RAUIICH 

WAVE  SPEED  WAVE  SPEEI 

(M/SEC)  (M/SIC ) 


e.eec 

0.000 

1527.061 

51.000 

51.000 

1523.539 

101 .003 

101.000 

1520.411 

152.000 

152.000 

1517.543 

202.000 

202.000 

1515.024 

253.000 

253.000 

1512.728 

333.ee? 

303. eee 

1510.725 

354.000 

354.000 

1508.915 

404 .000 

404 . 030 

1507.351 

455.000 

455.000 

1505.952 

505.000 

505.000 

1504.763 

556.000 

556.000 

1503.712 

606.000 

606. C00 

1502.636 

657.000 

657.000 

1502.085 

707.000 

707.000 

1501.477 

756.000 

758.000 

1500 .978 

606.0C0 

808.000 

1500.598 

859.000 

859.000 

1500.314 

909.000 

509.000 

1500.127 

960.000 

960.000 

1500.024 

1010.000 

1010.000 

1500.001 

1061.000 

1061.000 

1500.053 

1111.000 

1111.000 

1500.170 

1162.000 

1162.000 

1500.352 

1212. eec 

1212.000 

1500. 5e? 

1263.030 

1263.000 

1500.860 

1313.000 

1313.000 

1501 .216 

1364.000 

1364.000 

1501 .603 

1414.000 

1414. e00 

1502.024 

1465.000 

1465.000 

1502.492 

1515.000 

1515.000 

1502. 9S5 

1566.000 

1566.000 

1503.521 

1616.000 

1616.000 

1504.076 

1667.000 

1667. eee 

1504.669 

1717.000 

1717.000 

1505.276 

1768.000 

1768.000 

1505.919 

1818.000 

1818.00? 

1506.570 

1869.000 

iee9.000 

1507.254 

1919.000 

1919.000 

1507.943 

1970.000 

1970. tee 

1508.662 

2020.000 

2?20 . eee 

1509.383 

2071 .000 

2071.000 

1510.133 

2121 .000 

2121 .eee 

1510.860 

2172.000 

2172.000 

1511.655 

2222.000 

2222.000 

1512.426 

2273.000 

2273.000 

1513.223 

2323.000 

2323.000 

1514.013 

TOTAL  LOCAL  COMPRESS IOMAL 

DEPTH  DEPTH  WAVE  SPEED 

(M)  (M)  (M/SEC) 


2374.000  2374.000  1514. 828 
2424.000  2424.000  1515.635 
2475.000  2475.000  1516.465 
2525.000  2525. e00  1517.286 
2576.000  2576.000  1518.130 
2626.000  2626.000  1518.963 
2677.000  2677.000  1519.818 
2727.000  2727.000  1520.661 
2778.000  2778.000  1521.525 
28 28.000  2828. 0e0  1522.377 
2879.000  2879.000  1523.249 
2929.000  2929.000  1524.108 
2980.000  2980.000  1524.988 
3030.000  3030.000  1525.853 
3081.000  3081.000  1526.738 
3131.000  3131.000  1527.608 
3182.000  3182.000  1528.499 
3232.000  3232.000  1529.374 
3283.000  3283.000  1530.268 
3333.000  3333.000  1531.147 
3384.000  3384.000  1532.045 
3434. 0e0  3434.000  1532.926 
3485.000  3485.000  1533.827 
3535.000  3535.000  1534.712 
3586. eee  3586.  eee  1535. e'i6 
3636.000  3636.000  1536.503 
3687.000  3687.000  1537.408 

3737. eee  3737. eee  1538.297 
376e.000  3788.000  1539.205 
3838.000  3838.000  1540.095 
38&9.eee  3939. eee  1541.005 
3939.000  3939.000  1541.697 
3990.000  3990. eee  1542.807 
4040. eee  4040.000  1543.700 
4091.000  4091.000  1544.612 
4141.000  4141.000  1545.506 
4192 .000  4192. eee  1546.419 

4242. eee  4242. eee  1547.314 

4293.000  4293.000  1548.227 
4343.000  4343.000  1549.123 
4394.000  4394.000  155e.03? 
4444.000  4444.000  1550.934 
4495.000  4495.000  1551.349 
4545.000  4545.000  1552.746 
4596.000  4596.000  1553.661 
4646.000  4646.000  1554.558 
5000.000  5000.000  1560.917 


«. 

5eee.ee0 

5100.000 

e.eee 

lee.eee 

1550 .000 
1650.000 

1.540 

39 

3 

5iee.e0B 

5200.000 

e.eee 

lee.eee 

i65e .eee 
1680.000 

1.750 

38 

}  .  t  - 

5200. eee 

e.eee 

1703.000 

e.eee 

e.eee 

2.000 
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MODE  CHARACTERISTICS 


MODE 

NUMBER 

HATE  NUMBER 
(1/M) 

PHASE  VELOCITY 
(M/SEC) 

CROUP  TELOCITY 
(M/SEC) 

ATTENUATION  TERMS 

1 

.209300849 

1500.4991 

1500.0009 

.1000166031*01 

0. 

0. 

0. 

0. 

2 

.209231409 

1501.4919 

1500.0084 

.1000494261*21 

0. 

0. 

0. 

0. 

3 

.209094015 

1502.4785 

1500.0227 

.100081806E+01 

0. 

0. 

0. 

0. 

4 

.208957603 

1503.4594 

1500.0443 

.1001137332*01 

0. 

0. 

0. 

e. 

5 

.208822191 

1504.4343 

1500.0725 

.1001452231*01 

0. 

0. 

0. 

0. 

6 

.208687766 

1505.4034 

1500.1074 

•lC0176282E*ei 

0. 

0. 

0. 

0. 

? 

.208554294 

1506.3668 

1500.1480 

.100206944E*01 

0. 

0. 

0. 

0. 

a 


208421764 


1507.3247 


1500.1040 


1002371052*01 


.206200160 


1508.2770  1500.2476 


10  .208150474  1 500.2240  1500.3062 


11  .208020680  1510.1656  1500.3600 


0. 

0. 

0. 

0. 


.1002670542*01 

0. 

0. 

0. 

0. 


.100206 520 2*61 

0. 

0. 

0. 

0. 


.1003256532*01 


0. 

0. 


12  .207000757  1511.1021  1500.4385 


.1003543782*01 

0. 

0. 

0. 

0. 


13 

.207772716 

1512.0333 

1500.5118 

.1003827852*61 

0. 

0. 

0. 

0. 

14 

.207645512 

1512.9596 

1500.5901 

.1004108432*01 

0. 

0. 

0. 

0. 

15 

.20751P141 

1513.6809 

1500.6715 

.1004386182*01 

0. 

0. 

0. 

0. 

16 

.207393575 

1514.7975 

1500.7568 

.1004660942*01 

0. 

0. 
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vt"' 


•207268802 


.207144803 


.207021525 


206808936 


206776963 


206655527 


206534551 


206413938 


1515.7094 


1500 . 6467 


1516.6167 


1500.9353 


1517.5198  1501.0229 


1518.4190 


1501.1075 


1510.3146 


1501.1826 


1520.2074 


1501.2477 


1521.0979 


1501.2999 


1521.9867 


1501.3375 


III-ll 


.100493244E~ei 

0. 

0. 

0. 

0. 


.100 5202711*01 

0. 

0. 

0. 

0. 


.100  54  71 80  l+0'l 

0. 

0. 


.1C0574842E+01 


.1006010971+01 

0. 

0. 

0. 

0. 


.1006283771+01 

0. 

0. 

0. 


•10065599SE+01 


.1006840461+01 


■-i  >•<  * 


25 

.206293505 

1522.8745 

1501.3599 

. 10071 25561+01 
0. 

0. 

0. 

0. 

26 

.206173442 

1523.7620 

1501.3723 

.1007413761+01 

0. 

0. 

0. 

0. 

27 

.206053429 

1524.6495 

15ei.3?80 

.1007704131+01 

0. 

0. 

0. 

0. 

28 

.205933512 

1525.5374 

1501.3803 

.1007995611+01 

0. 

0. 

0. 

0. 

28 

.205813686 

1526.4255 

1501.3862 

.1008265842+01 

0. 

0. 

0. 

0. 

30 

.205683961 

1527.3140 

1501.3962 

.1008574651+01 

0. 

0. 

0. 

0. 

31 

.205574351 

1528.2026 

1501.4155 

.1008860271+01 

0. 

0. 

0. 

0. 

32 

.205454880 

1529.0812 

1501.4422 

.1009143241+01 

0. 

0. 

0. 

0. 

33 

.205335596 

1529.9796 

1501.4777 

.1009423061+01 

0. 

0. 

0. 
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229216505 


1530.6674 


1501.5261 


1009698211*2; 


205097665 


204979097 


224660623 


204742872 


204625265 


.204506017 


.204391136 


.204274661 


1531.7545 


1532.6405 


1533.5254 


1534.4086 


1535.2907 


1536.1709 


1537.0493 


1537.9256 


1501.5667 


1501.6534 


1501.7329 


1501.8176 


1501.9117 


1502.0126 


1502.1210 


1502.2427 


0. 

0. 

0. 

0. 


.100996862  £*  0fg] 

i: 


.101223695i.*0l| 


.10105008C£*0] 

0. 

0. 

0. 

0. 


.1010761141*01 

0. 

0. 

0. 

0. 


.1011018211*0] 

0. 

0. 

0. 

0. 


. 101127232£*0i 

0. 

0. 

e. 

0. 


. 101152 320E*01 

0. 

0. 

0. 

0. 


.1011768771*01 
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43 

.204158596 

1538.8001 

1502.3649 

.12120133SE+01 

0. 

0. 

0. 

0. 

44 

.204042940 

1539.6723 

1502.4916 

.10122  5560  £*-01 
0. 

0. 

0. 

0. 

45 

.203927701 

1540.5424 

1502.6258 

.1012494481+01 

0. 

0. 

0. 

0. 

46 

.203812890 

1541.4102 

1502.7650 

.101 2730781+01 
0. 

0. 

0. 

0. 

47 

.203698514 

1542.2757 

1502. 9ie4 

.Itl2964iel+01 

0. 

0. 

0. 

0. 

48 

.203584580 

1543.1388 

1503.0596 

.1213195252+01 

0. 

0. 

0. 

0. 

49 

.203471085 

1543.9996 

1503.2120 

.1013424361+01 

0. 

e. 

0. 

0. 

be 


203358235 


1544.8579 


1503.3699 


1C1365072E+01 

6711712771-18 


51 

.223245433 

1545.7138 

1503.5302 

.1013875351*01 
.  1763133562-16 

0. 

0. 

0. 

52 

.203133270 

1546.5673 

1503.6004 

.1014099081*01 
.278639986 E~1 5 

0. 

0. 

0. 

S3 

.203021547 

1547.4163 

1503.8591 

.1014319032*01  : 

.  38594 2629 E- 14 
.1418686222-19  . 
0. 

0. 

54 

.202010280 

1548.2666 

1504.0322 

.1£14S3655E*01 

.4665101972-13 

.3728873062-18 

0. 

0. 

55 

.202790455 

1549.1130 

1504.2025 

.1014754092*01 

.5597663062-12 

.6267555452-17 

0. 

0. 

56 

.202669068 

1549.9566 

1504.3776 

.1014969052*01 

.5869516652-11 

.7335975342-16 

0. 

0. 

57 

.202579124 

1550.7978 

1504.5567 

.1015181772*01 
.5596171622-10 
.7568131342-15 
.209550 5462-21 

0. 

56 

.202469616 

1551.6366 

1504.7334 

.1015394382*01 
.4639213452-09 
.7044867932-14  j 
.2067012932-20  j 
0.  ] 
3 

56 

.202360536 

1552.4730 

1504.9159 

I 

.1015604082*01  j 
.3786642372-08 

*  i 
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.5634450702-13 

.1845000261-19 

e. 


60 

.202251899 

1553.3069 

1505.1027 

.1815811411*01 

.2671648281-07 

.4507700911-12 

.1484906931-18 

0. 

61 

.202143694 

1554.1383 

1505.2897 

.101 601 7662-01 
.1692461881-06 
.3074457252-11 
.1873046962-17 
0. 

62 

.202035919 

1554.9674 

1505.4796 

.1016221332-01 

.9560404552-06 

.1873839932-10 

.6928957692-17 

0. 

63 

.201926572 

1555.7940 

1505.6708 

.1016420632*01 

.4815709192-05 

.1014223712-09 

.3973142862-16 

0. 

64 

.201821650 

1556.6183 

1505.8630 

.1016606032*01 

.2132295642-04 

.4835852142-09 

.2006870312-15 

0. 

65 

.201715140 

1557.4402 

1506.0469 

.1016748952*01 

.8224083202-04 

.2008611462-08 

.8830221702-15 

0. 

66 

.201608995 

1558.2602 

1506.1913 

.1016775702*01 

.2719492182-03 

.7153591902-08 

.3331376462-14 

0. 

67 

.201503035 

1559.0796 

1506.1797 

.1016562322*01 

.7543620001-03 
.2137689982-07 
.  1 054643281-1 3 

0. 


i 


III-16 


* 


68 

.201396738 

1559.9024 

1505.7987 

. 1 21 599440 E*01 
.171755595E-02 
.5246458421-07 
.2^431 53881-13 
0. 

69 

.201289147 

1560.7362 

1504.8423 

.101 51230*1*01 
.3183124311-02 
.1049447141-06 
.  582048921 E- 13 
0. 

70 

.201179063 

1561.5903 

1503.2999 

.1014219791*01 
.  4886953 58E-02 
.  1 742761 91 E-06 
.102691 550 i-1 2 
0. 

71 

.201065535 

1562.4720 

1501.4005 

.101 356461 E*£l 
.6474685361-02 
.250474437 E-06 
.1 57136109E-12 
0. 

72 

.200948140 

1563.3848 

1499.3953 

.1013Z6804E+01 
.775188068E-02 
.3263633371-06 
.  2 18504470 E- 12 
0. 

73 

.200826855 

1564.3290 

14 S7. 4303 

.1013292891*01 
.670599480 E-02 
.4002183821-06 
.2866468161-12 
0. 

74 

.200701852 

1565. 3e33 

1465.5573 

.101355180E*ei 
.9406033811-02 
.473667596 E-06 
•363790368E-12 
0. 

75 

.200573347 

1566.3061 

1493.7778 

.1013967041*01 

.99264877XE-02 

.5494125341-06 

.4534768431-12 

0. 

76 

.200441542 

1567.3361 

1492.0796 

.1014481181*01 
.103342803E-01 
.830322586 E-06 

Ill— 17 


.5603401571-12 

0. 


77 

.200306609 

1568.3916 

1490.4454 

.101 50  57471*01 
.1066548461-01 
.7 19147400 E-06 
.6900029761-12 
0. 

78 

.200168687 

1569.4726 

14  ee.  86i  5 

.10156?1351*ei 

.1095101291-01 

.8186767301-06 

.8495320261-12 

0. 

79 

.200027889 

1570.5773 

1487.3167 

.1016306351*01 

.1121080711-01 

.9318982701-06 

.1047959621-11 

0. 

80 

.196884311 

1571.7055 

1465.8021 

.1016951381*01 
.1145875721-01 
.10621 55771-05 
.1296960951-11 

0. 

81 

.196738030 

1572.8565 

1464.3112 

.1017598501*01 

.1170490211-01 

.1213336491-05 

.1611982061-11 

0. 

82 

.199589112 

1574.0301 

1482.8381 

.1018242231*01 

.1195654921-01 

.1390042331-05 

.2013221781-11 

0. 

83 

.199437612 

1575.2258 

1481.3786 

.1018877811*01 

.1221976191-01 

.1597884611-05 

.2527831021-11 

0. 

84 

.169283579 

1576.4433 

1479.9301 

.1019500651*01 

jf  !  1249996361-01 

.1843791061-05 

.3192343751-11 

l  •• 

I 


■V  a  .wt/f  :A 
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85 

.199127057 

1577.6825 

1478.4895 

i 

.i020ie?25E*ei j 
.1280156811-01 1 
. 21362 6972E-05! 
.4056106971-11 
*• 

86 

.198968081 

1578.9430 

1477.0548 

.102069410E*01 ' 
.1312S769eE-01 
.2485930451-05' 
.5186369551-11 

e. 

87 

.198806687 

1580.2248 

1475.6246 

.1021256971*01 
. 1 34 661 790 E- 01 
.29061 5464E-05 
.6675553731-11 
0. 

88 

.196642905 

1581.5277 

1474.1975 

.1021751621*01 
.1367624711-01 
.341379251 E-05 
•S65129359E-1 1 
0. 

89 

.198476764 

1582.8516 

1472.7729 

.102229369E+01 
.143100  51  5E-01 
.4030295321-05 
•112913784E-10 

e. 

90 

.198308291 

1584.1963 

1471.3504 

.1022757321+01 

.1478700111-01 

.47330202SL-05 

.1464498651-10 

0. 

91 

.198137513 

1505.5618 

1469.9297 

.1023176591*01 
•153151956E-01 
.57071 0999 E-05 
. 166641 127 E-10 
0. 

92 


.197964456 


93 


.19778914? 


i'lHil  i  ‘  I  >1 


1586.9476 


J. 


1588.3544 


111-19 


1468.5188 


1467.0843 


.1023544481*01 
.159015236E-01 • 
.6847930641-05 i 
.2624963151.-10 
0. 


.1023852951*01 j 
.1655379001-01 
.8264339751-05! 
.3531962971-10] 


■nour.-- 


1 


0 


04 

.107611612 

1580.7814 

1465.6809 

.1024002721*01 

.1728086711-01 

.1003307781-04 

.4791134631-10 

0. 

05 

.107431870 

1591.2287 

1464.2717 

.1024253261*01 

.1800275891-01 

.1225461271-04 

.6553394701-10 

0. 

06 

.197240079 

1502.6061 

1462.8680 

.1024322651*01 

.1500066761-01 

.1506097381-04 

.9039892641-10 

0. 

0? 

.197065945 

1594.1834 

1461.4717 

.1024287641*01 
.2001 69497 E-0 1 
.1662614821-04 
.1257687021-09 
0. 

96 

.196879813 

1595.6906 

1460.0851 

.l£24133S3E*ei 
.2115487161-01 
.2317972061-04 
. 176485002 E-e9 
0. 

90 

.196691623 

1597.2173 

1458.7107 

.1023846231+01 
. 224260 907E-01 
.2902449371-04 
.249769866E-09 
0. 

iee 

.106501416 

1598.7634 

1457.3514 

.1023410211*01 

.2384964721-01 

.365565745E-04 

.3564354241-09 

0. 

101 

.196300248 

1600.3264 

1456.0104 

.1022613331+01 

.2543021371-01 

.4630176461-04 

.512696418E-09 

0. 

102 

.196115166 

'  1601.9121 

1454.6005 

.1022048201+01 
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.2717526451-01 
. 56624 5566E-04 
.7426470161.-09 
0. 


103 

.165919230 

1603.5142 

1453.3940 

•102111696E+01 
.2908080361-01 
.7527507981-04 
.1083  1  2364  E-efi 
0. 

104 

.195721497 

1605.1342 

1452.1216 

.1020037591*01 
.3112740371-01 
.6636439651-04 
.  1 58'708820I-08 
0. 

105 

.195522024 

1606.7716 

1450.8713 

.1916851721*01 

.332727117E-01 

.1235035?eE-03 

.2332741731-08 

0. 

106 

.195320856 

1608.4266 

1449.6372 

.ltl763326E*01 
.35443 563ei-01 
. 15793 8989 E-03 
.343123654E-08 
0. 

107 

.195118033 

1610.0986 

1448.4073 

.1016494781*01 

.375304533E-01 

.2010113931-03 

.5036399191-08 

0. 

108 

.194913548 

1611.7877 

1447.1634 

.1015586311+01 

.3938899471-01 

.253760155E-03 

.7353556211-08 

0. 

109 

.194707370 

1613.4945 

1445.8814 

.  10150 8193E*ei 
.4085374141-01 
.31 677 7566 E-03 
. 1C6460663E-07 
0. 

110 

.194499417 

1615.2196 

1444.5342 

.1015150011*01 

.4170668431-01 

.3697837111-03 

.1523976651-07 

0. 
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111 

.194289566 

1616.9642 

1443.0975 

.1015911811*01 

.42ei60e43E-01 

.4716045961-03 

.2153030921-07 

0. 

112 

.194077657 

1618.7297 

1441 . 5553 

.1017403841*01 

.4157B6753E-01 

.5613645731-03 

.3000441481-07 

0. 

113 

.193863519 

1620.5177 

1439.9050 

.1019560331*01 

.4051771591-01 

.6571031891-03 

.412914595E-07 

0. 

114 

.193646991 

1622.3297 

1438.1573 

.1022227311*01 

.3696644821-01 

.7584807121-03 

.5626157251-07 

0. 

115 

. 193427945 

1624.1669 

1436.3331 

.1025200681*01 

.3718089911-01 

.8663619521-03 

.7616662301-07 

0. 

116 

.193296292 

1626. 03e2 

1434.4581 

.1026270271*01 

.3529941491-01 

.9833521711-03 

.ie2987eill-06 

0. 

117 

.192981992 

1627.9201 

1432.5580 

.1031252551*01 
. 3351 1 63231-0 1 
.1113949111-02 
.1396103941-06 
0. 

116 

.192755046 

1629.8366 

1430.6549 

.1034005371*01 

.3194482191-01 

.1264679161-02 

.1906309001-06 

0. 

ue 

.192525488 

1631.7801 

1428.7671 

.1036426091*01 

.3068567241-01 


.144444728E-02 

.2632382511-06 

0. 


120 

.102293378 

1633.7498 

1426. 9296 

.12384 42 83E*21 
.2879011621-01 
. 16652 5970E-02! 

.36681 «656E-26l 
0. 

121 

.192058790 

1635.7452 

1425.0957 

.1239989961*01 

.2929594991-ei 

.1943527351-22 

.5263698001-06 

2. 

122 

.191821854 

1637.7658 

1423.3391 

■  lk.412132  51*21 

.2923475661-01 
.23221 9925E-02 
.7667317751-06 


i 


i 


123 

.181582665 

1639.8105 

1421.655b 

.lf4143567E*01j 
•2P6418766E-G1) 
.277414720E-02! 
.11430  5629E-05 : 
0 . 

J 

124 

.191341384 

1641.8783 

1420.0675 

1 

.104U5321E*01 
.325642081 E-Kl 
.3427470 7ei-02; 
.1747967011-051 
0. 

125 

.191098197 

1643.9677 

1418.6031 

. 1242016321*01 
. 320 653672 E-01 
.427371211  E-02.1 
.2747105951-05 
0. 

126 

.192853340 

1646. e769 

1417.3035 

.123782595E*ei 
.3422587561-01 
.5480193711-02 
. 44433 1 928 E-05 
0. 

127 

.190607119 

1648.2032 

1416.2237 

.103430E21E*01 

•371319122E-01 

.7187229301-02] 

.739990724E-0» 


128 

.100339924 

1650.3435 

1415.4323 

.1029144311*01 

.4083620121-01 

.962723700E-02 

.1267161741-04 

e. 

12S 

.190112243 

1652.4936 

1414.9993 

.1022077671*01 
.4525962431-01 
.  131 10  5546E-0 1 
.2220644721-04 
0. 

130 

.160864627 

1654.6487 

1414.9534 

.1013164961*01 
.4996223721-01 
.17971 59611-01 
.3 £42 18 882 1-04 
0. 

131 

.189617555 

1656.8047 

1415.1902 

.1003302921*01 

.5401520041-01 

.2436660131-01 

.666*855341-04 

e. 

132 

.189371150 

1638.9605 

1415.3592 

.9647371371*00 

.5562635241-01 

.3168294921-01 

.1194990241-03 

0. 

133 

.189124817 

1661.1213 

1414. 6984 

.9905877801*00 

.5316828791-01 

.3927543581-01 

.1944081881-03 

0. 

134 

.188877208 

1663.2990 

1413.3935 

.9926649901*00 

.4652008041-01 

.4502084021-01 

.2969284601-03 

e. 

135 

.188626794 

1663.5071 

1420.9755 

.9695277061*00 

.3763? 7337 £-01 
. 4 8555 6247 £-0 1 
.4313653071-03 
0. 


136 


186372670 


1667.7540 


1408.2369 


.1 007458871*01 
.2911199631-01 
.5075986441-01 


6146014231-03 


,168114901  1670.0392  1403.6333 


0. 


.1C1284031E+01 
.2261149271-01 
.53279157*1-01 
.690223531  i-e  3 
0. 


187854397  1672.3551  1404.3040 


.1012985741*01 
. 187676860 E-01 
.5785801221-01 
.134679542E-02 
0. 


,187692781  1674.667*  1404.1781 


.1005445991*01 
.17884461 5E-01 
.662458968E-01 
.216131485E-02 
0. 


,187332432  1677.0148  1406.0389 


.987592312E*P0 
.2053689401-01 
. 7993181 66E-01 
.365567000 E-e2 
0. 


187076297  1670.3109  1409.9170 


. 95982 7914E*00 
.272574e54E-01 
.9758581121-01 
.6209867071-02 
0. 


186825894  1681.5617  1413.1647 


.93592 8275E+00 
.358656069E-01 
. 108970 632 E*00 
.0518025691-02 
0. 


186577179  1683.8033  1410.8728 


.940461261 E*00 
.3974096011-01 
. 9987731 52E-01 
. 1 182570021-01 
0. 


186321993  1686.1094  1403.2262 


.9714219861*00 

.3622317001-01 

.7593564861-01 

.1211497451-01 


III-25 


,ur,rj^i  I#  wv-.  .  *v 


e 


145 

.166056320 

1688.5170 

1395.2417 

.1003569881*01 

.3010335051-01 

.5436231551-01 

.116=410781-01 

0. 

146 

. 165781 564 

1691.6142 

1369.3634 

. 1024667361*01 
.2516534531-01 
.4il57e844E-01 
.lie697739i-01 
0. 

147 

.165500676 

1693.5748 

1385. 7072 

.1034436501*01 
.2215966551-01 
.3*6244351 1-01 
.  137*369731-01 
0. 

148 

. 1 85216616 

1696.1721 

1364.4528 

.1 03257264  E+01 
.2o8'»65405E-01 
.3351949571-01 
.1954800161-01 
0. 

149 

. 1P4 933632 

1698.7658 

1369.4109 

.1003922471*01 

.2110512771-01 

.3692474731-01 

.4330241471-01 

0. 


C.  DATA  SET  (2) 

The  following  printout  constitutes  the  data  set  for  the  second 
example  generated  by  NEMESIS. 


III-27 


Hl'KK  PKCFIIE,  Cl  -  1500,  Z1  -  1000,  EPS 
50. 


0074,  E  *  1233.33 


1.000 
0.002 
222.000 
404.000 
606.000 
628.006 
1010.000 
1212.000 
1414.000 
1616.000 
161S.000 
2C20.00C 
2222.000 
2424.000 
2626.000 
2829.000 
3030.000 
3232.00k 
3434.000 
3636.000 
3838 .000 
4040.000 
4242.000 
4444.000 
4646.000 


1527.081 
1515.024 
1507.351 
1502 .836 
1500.598 
1500.001 
I5e0.587 
1502.024 
1504.076 
1506.570 
150 P.3 83 
1512.426 
1515.635 
1518.963 
1522.377 
1525.853 
1529.374 
1532.926 
1536.503 
154£ .095 
1543.700 
154  *7. 314 
1550.934 
1554. 55P 


51.000 
253.000 
455.000 
657. eee 

ess.ep® 
1061.000 
1263.000 
1465.000 
1667.000 
1669.000 
2071.000 
2273.000 
24'»5.0e0 
2677.000 
2879.000 
3e61.000 
3263.000 
3485.000 
3687. £00 
3889.000 
4091.000 
4283.000 
4495.000 
5000.000 


1523.539 
1512.726 
1505.952 
1502.085 
1500.314 
1500.053 
1500.680 
1502.492 
1504.669 
1507.254 
1510.133 
1513.223 
1516.465 
1519.816 
1523.249 
1526  .73e 
1530.268 
1533.627 
1537 .408 
1541. zes 
1544.612 

1548.227 

1551 .849 
1560.917 


101.000 
303.ee? 
505.000 
7C7.000 
909.000 
1111 .000 
1313.000 
1515.000 
1717.000 
1816.000 
2121.000 
2323.040 
2525.000 
2727.040 
2920.000 
3131. eu' 
3333.000 
3535.000 
3737. ect 
3938.00k. 
4141.000 
4343.000 
4545.000 


1520.411 
1510.725 
1504  .760 
isei.477 
1500.127 
1500.170 
1501.216 
1502.985 
1505.276 
1507.943 
1510.860 
1514.013 
1517  .286 
1520.651 
1524. 10S 
1527  .50? 
1531.14? 
1534  .712 
1538.297 
1541.897 
1545.506 
1549.123 
1552.746 


152.000 
354 .000 
556.000 
758.000 
960.000 
1162.000 
1364. eee 
1566.000 
i76e.000 
1970.000 
217?. e00 
2374.000 
2576.000 
277t .000 
2960.000 
3162.000 
3364.000 
3586.000 
3768.000 
3990.000 
41 92.000 
4394. Z0Z 
4596.000 


1517.543 
1 508 .915 
1503.712 
1500.978 
1500.024 
1500.352 
1501.603 
1503.521 
1505.919 
1506.662 
1511.655 
1514.826 
1513.130 
1521.525 
1524.966 
1523.499 
1532.045 
1535.616 
1533.205 

1542.607 

1546.419 

1550.037 

1553.661 


1.500 

5020.000  1550.000  5100.000  165e.000 

1.750 

5i00.0?e  lfse.eeo  5200.000  ieee.000 

2.000 

5200.000  5000.000  5200.000  2700 .000 


I 


D.  PRINTED  OUTPUT  (2) 


The  following  printout  is  the  output  of  data  set  (2)  from  NEMESIS. 


*■ 
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MIMESIS  ARL-UT  NORMAL  MOLE  MOLEL 
VERSION  1.2 


monk  profile,  ci  *  isee.  zi  •  1200,  ips  «  .0074,  b  -  1233.33 


FREQUENCE  (HZ)  50.000 

NUMBER  Or  EIISTIN3  MOLES  292 

FIRST  MOLE  COMPUTED  1 

LAST  MODE  COMPUTED  26? 

DEPTH  (M)  AT  KHICH  MODES  ARE  MATCHED  1010.000 


NO  OPTIONS  SELECTED 


i 

* 

i 

t 


l 


) 


I 


} 

I 


GEOACOUSTIC  PROHLE 


i 


TOTAL 

LEPTH 

(M) 

0.020 

0.000 

51.000 

51.000 

101 .000 

121.000 

152.002 

152.000 

202.000 

202.000 

253.020 

253.200 

303.000 

303.000 

354.000 

354.000 

404.200 

404.000 

455.000 

455.000 

505.000 

505.000 

556.000 

556.000 

606.000 

606.200 

657.000 

657.000  ■ 

727.020 

727.000 

758.000 

758.000 

808.000 

806.000 

859.200 

859.200 

909.000 

909.020 

960.000 

960.000 

1210.002 

1010.000 

1061.000 

1061.200 

1111 .000 

111 1.000 

1162.000 

1162.200 

1212.000 

1212.000 

1263.000 

1263.000 

1313.002 

1313.202 

1364.000 

1364.000 

1414.000 

1414.000 

1465.000 

1465.222 

1515.000 

1515.000 

1566.000 

1566.000 

1616.020 

1616.202 

1667.000 

1667.200 

1717.000 

1717.000 

1768.000 

1768.220 

1818.000 

1P18.000 

1869.000 

1869.000 

1919. 02e 

1919. oeo 

1972.000 

1970.000 

2020.000 

2020.000 

2271 .022 

2071.200 

2121.000 

2121.000 

2172.000 

2172.000 

2222.000 

2222.200 

2273.000 

2273.200 

2323. 00e 

2323. e00 

COMPRESS IONAI  SHEAR 

VATE  SPEED  WAVE  SPEFD 

(M/SEC)  (M/SEC) 


1527.061 
1523.539 
1522 .411 
1517.543 
1515.224 
1512.723 
1512.725 
1508.915 

1507.351 
1505.952 
1504.760 
1523.712 
1502. e36 
1502.085 
1521 .477 
1500.978 
1500.598 
1502.314 
1500.127 
1500.024 
1580. eei 
1500.053 
1500.170 

1500.352 
1500.587 
1500.880 
1501.216 
1501.603 
1502.024 
1522.492 
1502.985 
1503.521 
1504.076 
1504.669 
1505.276 
1525.919 
1506.570 
1507.254 
1507.943 
1508.662 
1509. 3e3 
1510.133 
1510.860 
1511.655 
1512.426 
1513.223 
1514.013 


RATDEIGH  SAMPLING 

WAVE  SPEED  DENSITY  POINTS 

(M/SEC)  (C/CC ) 


1.000  1921 
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•  •  '  JC-  •  •>  A"  :  - 


2374.000 

2374.000 

1514.828 

2424.000 

2424. 000 

1515.635 

2475.000 

2475.000 

1516.465 

2525.000 

2525.000 

1517.286 

2576.000 

2576.000 

1518.130 

2626.000 

2626.000 

1518.963 

2677.000 

2677.000 

1519.818 

2727.000 

2727.000 

1520.661 

2778.000 

2778.000 

1521.525 

2828.000 

2628. eee 

1522.377 

2679.000 

28^9.000 

1523.249 

2929.000 

2929.000 

1524.108 

2960.000 

2960. eee 

1524.968 

3030.000 

3e30.600 

1525.853 

3081.000 

3061.000 

1526.738 

3131. e?0 

3131 .000 

1527. 6ee 

3162.000 

3182. eee 

1526.499 

3232.000 

3232.000 

1529.374 

3283.000 

3283.000 

1530.268 

3333.000 

3333.000 

1531.147 

3384.000 

3384.000 

1532.045 

3434. 0e? 

3434.600 

1532.926 

348  5.000 

3485.000 

1533. e27 

3535.000 

3535.000 

1534.712 

3566.000 

3586. C00 

1535.616 

3636.000 

3636. e00 

1536.503 

3687.000 

3687.000 

1537. 40B 

3737.000 

3737. eee 

1538.297 

3786.000 

3788. e00 

1539.205 

3838.000 

3838.000 

1540.095 

3869 .000 

3889. eee 

1541.005 

3939.000 

3939.000 

1541.897 

3990. 000 

3990.000 

1542.807 

4040.000 

4040. eee 

1543.700 

4091.000 

4e?i.000 

1544.612 

4141.000 

4141.000 

1545.506 

4192  .eee 

4192. eee 

1546.419 

4242.000 

*242.000 

1547.314 

4293.000 

4293.000 

1548.227 

4343.000 

4343. eee 

1549.123 

4394 .000 

4394. 000 

1550.037 

4444.000 

4444.000 

1550.934 

4495.000 

4495.000 

1551.849 

4545.000 

4545.000 

1552.746 

4596.000 

4596.000 

1553.661 

4646.000 

4646.000 

1554.558 

_  _  _  — 

5000.000 

5000.000 

1560 .917 

2 

5000.000 

0.000 

1550.000 

1.500 

30 

_  _  _  _ 

5100.000 

100.000 

1650.000 

3 

5100.000 

e.eee 

1650.000 

1.750 

39 

_  -  _  - 

52u0 .000 

100.000 

1680.000 

4 

5200.000 

0.000 

5000. 000 

2700  .eee 

2501 . 555 

2.00.0 

MODE  CHARACTERISTICS 


MODE 

NUMBER 

WAVE  NUMBER 

(1/M) 

PHASE  VELOCITT 
(M/SEC) 

CROUP  VELOCITT 
(M/SEC) 

attenuation  terms 

1 

.208369649 

1600.4991 

1500.0009 

.100£l 6603 1+01 

0. 

0. 

0. 

0. 

2 

.209231409 

1501.4919 

1500.0084 

.  10004 9429E+ 01 

0. 

0. 

0. 

0. 

3 

.209094015 

1502.4785 

1500.0227 

. 1 00081 806 E+01 

0. 

0. 

0. 

0. 

4 

.206857603 

1503.4594 

1500.0443 

.1001137331+01 

0. 

0. 

0. 

0. 

5 

.206622191 

1504.4343 

1500.0725 

.1001452231+01 

0. 

0. 

0. 

0. 

6 

.206687766 

1505.4034 

1500.1074 

.1001762821+01 

0. 

0. 

0. 

e. 

7 

.206554294 

1506.3668 

1500.148b 

.1002068441+01 

0. 
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.208421764 

1507.3247 

1500.1949 

.1102371952*01 

0. 

0. 

0. 

0. 

.209290160 

1508.2770 

1500.2476 

. 1002670 54E-01 
0. 

0. 

0. 

0. 

.208159474 

1509.2240 

1500.3062 

.  1 002 96520 E* 01 
0. 

0. 

0. 

0. 

.208029680 

•  1510.1656 

1500.3690 

.1003256531*01 

0. 

e. 

0. 

0. 

.207900757 

1511.1021 

1500.4385 

.10035437SE-01 

0. 

0. 

0. 

0. 

.207772716 

1512.0333 

1500.5118 

. 100382785E*01 

e. 

0. 

0. 

0. 

.207645512 

1512.9596 

1500.5901 

.  1004108431+0 1 
0. 

0. 

0. 

0. 

.207519141 

1513.8605 

1500.6715 

.  100438618E+01 
0. 

0. 

0. 

0. 

.207393575 

1514.7975 

1500.7568 

.100466094E+01 

0. 


z. 

e. 


17 

.207268802 

1515.7094 

1500.8467 

•100493244E+01 

e. 

e. 

0. 

e. 

ie 

.207144803 

1516.6167 

1500.9353 

.1005202711*0] 

0. 

e. 

e. 

0. 

19 

.207021526 

1517.5198 

1501.0229 

.1005471902*0] 
0  . 

0. 

0. 

0. 

20 

.206898936 

1518.4190 

1501 .la^s 

.  I  .,05740422*0! 
0  . 

0. 

0. 

e. 

21 

.206776963 

1519.3146 

1501 . 1826 

.1006010972*01 

0. 

0. 

e. 

0. 

22 

.206655527 

1520.2074 

1501.2477 

.10262  9377  2*0 
0 . 

0. 

0  . 

0. 

23 

.206534551 

1521 .0979 

1501 .2999 

.1006559952*0 

0. 

0. 

0. 

e. 

24 

.  206413938 

1521 .9867 

1501.3375 

.  1006c  4046  2*0 
0. 

0. 
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25 

.206293595 

1522.8745 

1501.3596 

.1007125561-01 

0. 

0. 

0. 

0. 

26 

.206173442 

1523.7620 

1501.3723 

.100741378E-01 

0. 

0. 

0. 

0. 

27 

.206053429 

1524.6495 

1501.3780 

.10077241 3  £4-01 
0. 

0. 

0. 

0. 

26 

.205933512 

1525.5374 

1501.3603 

.109799561 i+01 
0. 

0. 

0. 

0. 

29 

.205813668 

1526.4255 

1501.3662 

.10032S564E-01 

0. 

0. 

0. 

0. 

30 

.205693961 

1527.3140 

1501.3962 

•100857465E+01 

0. 

0. 

0. 

0. 

31 

.205574351 

1528.2026 

1501.4155 

.  1«;2F8602?I*01 
0. 

0. 

0. 

0. 

32 

.205454890 

1529.0912 

1501 .4422 

.100914324E»ei 

0. 

0. 

0. 

0. 

33 

.205335596 

1529.9796 

1501.4777 

.100942306£*01 

e. 

0. 

0. 
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20521(505 


205097665 


204979097 


204860623 


204742872 


.  204625265 


.204508017 


.204391138 


.204274661 


1530.8674 


1531.7545 


1532.6405 


1533.5254 


1534.4088 


1535.2907 


1536.1709 


1537.0493 


1537.9258 


1501.5261 


1501.5867 


1501.6534 


1501 .7309 


1501.8176 


15ei.9117 


1502.0126 


1502.1210 


1502.2427 


.1009698211*21 

0. 

0. 

0. 

0. 


.1009968621*01 

0. 

0. 

0. 

0. 


.1010236951*01 

0. 

0. 

0. 

0. 


.1010500891*01 

0. 

0. 

0. 


.1010761141*01 


.101 101 821 £*01 

0. 

0. 

e. 

e. 


.1011272321*01 

0. 

0. 

0. 

0. 


.1011523201*01 

0. 

0. 

0. 

0. 


.1011766771*01 


r 


43 


44 


45 


46 


47 


4e 


49 


I  ! 


k  50 


'  I 

s 


r  I 


.204158598 


.204042940 


.203927701 


.203812890 


.203698514 


.203584580 


.203471035 


.203358035 


1538.8001 


1539.6723 


1540.5424 


1541. 4iea 


1542.2757 


1543.1388 


1543.9996 


1544.8579 


1502.3649 


1502.4916 


1 5 *>2 . 6258 


1502.7650 


1502.9104 


1503.0596 


1503.2120 


1503.3699 


III-38 


0. 

e. 

e. 

0. 


.101201335E-01 


0. 

0. 

0. 

0. 


.1012255601+01 

0. 

0. 

0. 

0. 


.1012494481+01 

0. 

0. 

0. 

0. 


.101273078E-01 

0. 

0. 

0. 

0. 


.1012964101+01 


0. 

0. 

0. 

0. 


0. 

0. 

0. 


.1013195251+01 


0. 

0. 

0. 

0. 


.1013424361+01 


.1013650721+01 

.«7U7127?E-16 


0. 

0. 

0. 


~T 


.203245433 


.203133270 


.203021547 


.202910280 


.202799455 


.202689066 


.202579124 


.202469616 


1545.7138  1503.5302 


.1013875352*01 

.1763133582-16 


0. 

0. 

e. 


1646  5673  1503.6904 


.1*14099062*01 

.2786399862-15 

0. 

0. 

0. 


1547.4183  1503.8591 


.It 1431903E*01 
.3859426291-14 
.1416686221-19 
0. 

0. 


1548.2663 


1504.0322 


1549.1130 


1504.2025 


. 1 01 4536552*01 
.4365201972-13 
.3726673062-18 
0. 

0. 


.1014754092*01 

.5597693061-12 

.628755545E-17 

0. 

0. 


1549.9566 


1504.3776 


,1» 14969052*01 
.5669516952-11 
.7335975341-16 
0. 

0. 


1550.7978 


1504.5567 


.1015181772*01 
.5596171222-10 
.7566748511-15 
.46597 5273E-20 
.5233411692-20 


1551 .6366 


1504.7334 


.1015394382+01: 
.4o392134ci-09 
. 704550 649E-14 
.4711025062-19 
. 53442 9448E-19! 


202360538 


1552.4730 


1504.91 59 


.101 56040ei-0l! 
.3706642372-08! 


.5935011931-13 
.43121 6473E-18 
.4694330 571-18 


.202231899  1553.3069 


.202143694  1554.1383 


.202035919  1554.9674 


.221P28572  1555.7940 


.201621650  1556.6163 


.201715140  1557.4402 


.201606995  1558.2602 


.201503035  1559.0796 


1305.1027  .1615811411*01 

.2671648281-07 
.45081 60201-12 
.3560857851-17 
.4043664961-17 


1505.2897  .1016017661*01 

.1692461881-06 

.3074794861-11 

.2641641701-16 

.3001349701-16 


1505.4796  .1016221331*01 

.9580404581-06 

.1874061811-10 

.1752165041-15 

,19S17nlZZt-15 


1505.6708  .1016420631*01 

•481570921E-05 
•  U 143  5324  E-09 
.1032652561-14 
.117446481 E-14 


1505.8630  .1 £1660 6031*01 

.2132295651-04 
.483651 e571-0S 
.5364431561-14 
.6104223581-14 


1506.0469  .1016748951*01 

.8224083251-34 
.202  691028E-te 
.2429089851-13 
.2765482141-13 


1506.1913  .l£167757tl*Cl 

.2719492231-03 

.7154741301-08 

.9437638231-13 

.1075005951-12 


1506.1797  . 101656232 1’*01 

.7543620021-03 

.2138361171-07 

.307P27079E-12 

.3509261831-12 


201396736 


1559.9024 


1505.7967 


.1015994401+01 
.1717555961-02 
.5247444021-07 
.6262517261-12 
.9421 10374 £-12 


.201289147 

1560.7362 

1504.8423 

.1C1512304E+01 
.3183124351-02 
.1248660621-66 
. 16109 564 5E- 11 
.2065960941-11 

.201179063 

1561.5903 

1503.2999 

.1014219792-01 
.4F6695365E-02 
. 17431 4758E-06 
.3306246902-11 
.3773807641-11  1 

.201065535 

1562.4720 

1501.4005 

.101356461 £+01 
.6474685472-02 
.2505349031-06 
.5246662591-11 
.5991916842-11 

.200948140 

1563.3848 

1499.3953 

. 101 32 6804 E+ 01 
.7751860832-02 
.3264496352-06 
.7585159411-11 
.6657506512-11 

.200626855 

1564.3290 

1497.4303 

.1013292892+01  . 
.6705994992-02  : 
.4003346022-06 
.103,’37566i-10 
.1186099752-10  j 
3 

.200701852 

1565.3033 

1495.5573 

.1013551601+01  ! 
.9406034072-02  j 
.47381 8841I-e6  j 
.1376469602-10 
.1574765182-10 

.200573347 

1566.3061 

1493.7778 

.1013967042+01 
.9928488002-02 
.5496065212-06 , 
.1799541062-10 
.2060075462-10 

.200441542 

1567.3361 

1492.0796 

.1014481162+01 

.1033428062-01 

.6305732632-06 

12743246a 1-27 


94 

.197611612 

1589.7814 

1465.6810 

.12240  92492*01 
.1728094632-01 
.1015247362-04 
.20239534^2-07 
.2351777062-07 

95 

.197431879 

1591.2227 

1464.2717 

.  10242 5265 E*ei 
.1809288892-01 
.1247736232-04 
. 404681550 E-07 
.4706680452-27 

96 

.197249979 

1592.6961 

1462.8661 

.1224321842*01 

.1900089452-01 

.1552579312-04 

.912920042E"07 

.U627P6B3E-06 

97 

. 197065946 

1594.1834 

1461.4720 

.1024285712-01 
.2001738552-01 
. 16802 7994 r-04 
. 25264 9992 E-06 
•294408996E-06 

96 

.196879814 

1595.6906 

1460.0861 

.102412699E+01 
.2U558e96E-ei 
.2773960012-04 
.1085553742-05 
.  126621 723 E-0 5 

99 

.196691630 

1597.2173 

1458.7316 

.1023685141*01 

.2243210252-01 

.150844191E-03 

.3275027332-04 

.3323970012-04 

100 

.196640695 

1597.6310 

1599.5306 

.2069164682-03 

.262987246E-04 

.COfcJOrtijji  WT 

.be90143732*£0 
.2255932472 *00 
. 2634700 87 E*00 
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.2542788491-01 
.502661359E-04 
.1460581621-05 
.173221 67HE-05 


103 

.196115164 

1601.9122 

1454.6904 

.1022247621*21 
.2717326191-31 
.6066233131-04 
.8228410891-06 
.963686894 £-06 

104 

.195919228 

1603.5142 

1453.3937 

.lf2111765E*01 

.2927879111-01 

.’’620289821-34 

.6131797701-26 

.7186928071-06 

105 

.195721494 

1605.1342 

1452.1212 

.1020039011*01 
.3112522351-01 
.9666215731-04 
.5325517241-06 
.625029540 E-06 

106 

.195522021 

1606.7718 

1450.8708 

.1018853851*01 

.3327025611-21 

.1236073901-03 

.5078877881-06 

.5967260121-06 

107 

.195320856 

1608.4266 

1449.6365 

.1017636211*01 
.3544074741-01 
.1576751861-03 
.5151219391-06 
.6256670 49 E-26 

108 

.195118029 

1610.0986 

1448.4064 

.1016496601*01 
.3752721711-01 
.20031 7095E-23 
.5445614271-06 
.6412226931-06 

109 

.194913543 

1611.7876 

1447.1623 

.  101 55911 8E+ei 
.3938531251-01 
.2525447341-03 
.5915865771-06 
.6973755631-06 

110 

.194707363 

1613.4945 

1445. e799 

. 12150 8799E*01 
.4064963571-01 

.3148455661-83 

.653399394E-06 

.7711160001-06 


III-45 


.•v^An  *  i 


111 

.194499409 

1615.2196 

1444.5324 

.K1515739E+01 
.4176253061-01 
.3869557031-03 
.72774  e766i-k'6 
.6596475661-06 

112 

.194289556 

1616.9642 

1443.0953 

.101 5920651*01 
.420133996E-01 
.4676373261-03 
.8125916651-06 
.961 20 9747 E-k 6 

113 

.194077645 

1618.7298 

1441.5526 

.101  '’414331+01 
.41573S804F-01 
.555968427 1-0 3 
.90633621 3 t-06 
.107336687E-05 

114 

.193863504 

1620. 51 ?e 

1439.9018 

.  1 01 £572771*01 
. 40512 6663E-01 
.«4S932392E-03 
.3008436761-05 
. 1 19571665E-05 

115 

.193646974 

1622.3296 

1438.1534 

.1022242191*01 

.3896156221-01 

.7*S110615E-03 

.112«303918£-05 

.132966417E-05 

116 

.193427923 

1624.1671 

1436.3284 

.1025216751+01 
.3717579771-01 
.65424  7769  E-k  3 

.1244471291-^5 

.1479215241-05 

117 

.193206266 

1626.0304 

1434.4525 

.102829266E+01 
.352 9377461-01 
. 967 763097 E-03 
.1 36754642 E-0 5 
.1651356071-05 
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.1327956242-01 

.200514735E-01 

.1730087951-31 

.1084553472-02 

.2^14346792-02 

234 

.149950338 

2094.9630 

1111. e2l5 

.1334261342-01 

.1962116662-01 

.1676343482-01 

.1056844532-02 

.1966462232-02 

235 

.149381466 

2103. e672 

1107.4944 

.1339727911-01 
.1959523492-31 
. 166277073 E-01 
.1056331382-02 
.2310310072-02 

236 

.14679ee30 

2111.3020 

1103.2378 

.1344379702-0 1 
.1995653472-31 
.1689168212-0 1 
.1083063222-02 
.2381701332-32 

237 

.146211769 

2119.6648 

1099.0567 

.13431 91052-01 
.2071415612-31 
.1756438332-01 
.  1 134642S5l-£  2 
.2207467472-02 

236 

.147620654 

2128.1525 

1094.9626 

.1351093522-01 

.2166156472-01 

.1868319222-01 

.1216079932-02 

.2395759662-02 

239 

.147025903 

2136.7613 

1090.9722 

.1352968442-01 

.2349114722-01 

III-60 


.4 i 32 67064 1-3 2 


144620962  2172.2941  1076.3583  .  1348 82386 1*0 1 

.3435496811-01 
.326  14  55211-.;  1 
.2229336351-32 
.47345cr>0''l-02 


144015772  2181.4226  1072.9455  . 13461 3971 1 -5 1 

.375t73422i-0l 
.  3r~4326?2 1 -3 1 
.2541244751-32 
.0445912241-32 


143409211 


2190.6491 


1065.4089 


. 134  4  61 61 61  -  0  1 
.  396r'?3r'64  1-0  1 
.4043210201-31 
.2632152791-32 
.6161712721-02 


142600351  2199.5894  1065.5397  . 1 34654934 t-2 1 

.4055627241-01 
. 428302^301-01 
. 3.  4 652 604 1-6 2 
.6734677381-02 


142107536 


2209.4712 


1061 .1708 


•1354040431*01 
.35664252 61-3 1 
.4342607591-31 
.3138117921-02 
.7,  531 67151-3  2 


III-61 


248 

.141568761 

2219.1285 

1056.2721 

.1362772e5I+31 
.3717914371-81 
.4212287451*01 
.3296334841-02 
.7*  88 1 4 634 E “82 

249 

.142942278 

2228.9925 

1050.9633 

.1375976911*81 
.  337391 750 1-01 
.39421 6716E-21 
.2951472051-32 
.6  87134462  E"82 

2  50 

.142306988 

2239.8851 

1845.4336 

.139343701 1*8  1 
.3338997761-01 
.368 "34193E-1 1 
.2754400791-02 
.6534356691-02 

251 

.139662538 

2249.4166 

1339.8502 

.1404729971*01 
.2676696031-31 
.3..  "2  91 739 1-81 
.2551631681-32 
.6174387021-32 

252 

.135829124 

2259.9982 

1834.3281 

.141 58 1391 E*0 1 
.2411795941-01 
. 297921 9661-31 
.2374622061-22 
.5o6571237i-02 

253 

.138347209 

2278.8831 

1828.8963 

.1429929331*31 
.2217326261-81 
.27446496''  1-81 
.2236657641-02 
.56511 91341-32 

254 

.137677379 

2281 .8513 

1023.6333 

.1440386761*31 
. 2 1  943  92  6c 1-0 1 
.257,424241-01 
.2153503121-82 
.556371 6501-01 

255 

.137032252 

2293.1251 

1318. 4434 

.1449435661-01 
.2039262031-01 
.2466369241-.  1 
.2112656761-02 
.6587974341-0  2 

2  56 

.  136316448 

2324.6321 

1313.4263 

.1456992651-81 
.2 04 93 678 cl- 31 

2127001791-02 

57* a9e216E-ir 


257 

.135626616 

2316.3541 

1008. 5741 

.1463106611*01 
.2 124 5709? £-01 
•244250358E-01 
.2196712391-02 
.611°’ 08361-02 

258 

.134931451 

2328.2879 

ie03.9107 

. 1467*301 51*01 
.2269399661-01 
.2525193261-01 
.2327?  57671-02 
.6662040601-02 

25£ 

.134231737 

2340.4246 

999.4736 

.147?37-7^E*.'l 
.2439360691-01 
. 2  676 1 99661-0 1 
.2525125501-02 
.74425691 51-02 

262 

.133528373 

2352. 752S 

995.3066 

.14r,1149231-0i 

.2791700221-01 

.2699299491-01 
,£”97761441-0? 
.Prw  231312 1-02 

261 

.132822356 

2365.2589 

991 .43to7 

.1469672921-01 
.  3 177r,P560  £-0 1 

.3132873051-01 

.3146090361-02 

.9^73365501-0? 

262 

.13211466? 

2377.9269 

967.8463 

.1466693791- 21 
.  3631661c! !-.  1 
.353?c5P9el-01 
.  3  554  539931-02 
.  1 1 53c  96591-0 1 

263 

.  131425342 

2390.7539 

964 . 3653 

.  1463437?!  i*,.  i 
•  412303  4  7 1 “ . 1 
.  V-9392609 1-01 
.  59”514790i-02 
.1337003421-01 

264 

. 13? 696619 

2403.7401 

960.7390 

.It 6 16 3 591 :  - 0 1 
.4t9  6 ”116  7 1-01 
.41”c640T:'i-tl 
.431,-573631-0? 
.1506  3063  6 1-01 

II 1-63 
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265 

.129963376 

2416.9188 

976.4738 

.1465?  50602-01 
.4695432971.-01 
.427044372^-01 

.44^6330011-02 

.1626631501-01 

266 

.129265173 

2430.3473 

971.2652 

.1*75020062-01 

.4625591751-01 
.4159994661-01 
.4362419361-22 
.  1 661 6’7381  2-01 

26? 

.128536102 

2444.8945 

965.1266 

.1451219112-01 
.4313397181-01 
.3845946181-01 
.4  0644  S31 1 2-0  2 
.1612921011-01 

266 

.127799573 

2458.2184 

95S.3822 

.1511060612-01 

.3865136612-01 

.3433073112-01 

.3623556962-02 

.1510546962-01 

269 

.127048349 

2472.7536 

951 .4273 

.1  531'748332*0l 
.3395513472-01 
.3, 191692*1-01 
.3166657851-02 
.1352745972-01 

270 

.126284389 

2467.7126 

P44.5477 

.1551300431-01 
.  28’7970730  2-0 1 
.2665516932-01 
.2762297632-02 
.12981 6S9« 1-01 

271 

. 12E5?e352 

2523.0945 

937.8912 

.1566656131-01 
.26491 ’’2732-01 
.2556721202“. 1 
.2438305242-02 
.1212643372-01 

272 

.124721191 

2516.8524 

931.5203 

.15541 57652-01 
.2409768361-01 
,221653243t-.‘  1 
.2188139262-02 
.117.3654651-01 

273 

. 123923950 

2535.0973 

925.4706 

.1557157552-01 

.  225'r51 6451-01 
.2123973552-01 
.2035762762-02 


III-64 


e»s  *v  <%  rj  cw  pw  rw  ts>  r-j  c*»  e*  r-  ®  c«  n  w  c«  ^  ^ 


274 


.11761  e?79E.-ei 


.123117726  2551.6961 


919.7919 


275 


. 122333747 


2568.6806 


914.5814 


276 


.121483529 


2586.0235 


912.2222 


277 


.120659138 


2603.6923 


906.4239 


278 


.119833554 


2621.6302 


904.3056 


279 


.119011027 


2639.7492 


904.3796 


.116196866 


2657.9323 


907.2595 


261 


.117395322 


2676.0700 


912.9244 


.116608756 


2694. 13K 


921.0382 


. 162774 626E-01 
.2166267451-01 
.2il?S?e6cE-0’ 
.1944139661-32 
.1227016301-01 


.161558S53E-0* 

.2203107922-01 

.2606464922-01 

.1919001262-02 

.1336470442-01 


.161 991 6981-01 
.231 160474E-0 1 
.2408356561-01 
.  1 95990  624 1-02 

. 1 53291 3652-0 1 


.  1  cl  93  97071-, 0* 
.26^436  6 04  2-01 
.275^6 53151-01 
.2066655921-02 
.  1  o4  92  i  -0  1 


.16119 30 0rl* 01 
.69004000.  1-cl 
.33029370  52-0 1 
.2243409811-02 
.23517510"1-L i 


.159503r'55i»01 
. 344046784 1-01 
.4k c4036o6i-,  1 
.24547r'186i-02 
.3133335441-01 


. 156ol 61 161*0 1 
.4099526681-01 
•  5i 34  6360  71-0 1 
.2596834341-02 
•42?030”02i-01 


.  1 53r>0  765r,l  •  0 1 
.4615794341-01 
.5777466831-01 
.2435262672-0  2 
.5697238291-01 


.1486709921-01 


.441'7926991-01 
.5- 9k 347981-01 
.1672899642-02 
.9580467401-01 
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PHYSICAL  DESCRIPTION 


PLMODE  IS  A  COMPUTER  PROGRAM  DESIGNED  FOR  USE  WITH  THE  ARL :UT  NORMAL 
MODE  MOOEL*  NEMESIS.  THIS  MODEL  ASSUMES  A  HORIZONTALLY  STRATIFIED 

deep  ocean  with  single  channel  profiles  ano  multiple  fluid  sediment 
layers  overlying  a  substrate,  sound  speed*  in  the  water  and 

SEDIMENT  LAYERS  VARY  WITH  DEPTH.  THE  DENSITY  IS  CONSTANT  WITHIN 
EACH  LAYER.  THE  LAST  LAYER  IS  A  HOMOGENEOUS*  SEMl-lNFiNlTE 
FLUID  OR  SOLID  SUBSTRATE*  i.E.*ThE  COMPRESSlONAL  AND  SHEAR 
WAVE  SPEEDS  AND  OENSITY  REMAIN  CONSTANT  WITH  DEPTH.  SOUND  SPEEDS 
ANO  DENSITIES  CAN  BE  DISCONTINUOUS  AT  LAYER  INTERFACES.  THE 
ATTENUATIONS  IN  THIS  ENVIRONMENT  ARE  DEFlNFD  IN  PLMODE.  ATTENUATION 
IN  THE  WATER  ANO  THE  SHEAR  AND  COMPRESSlONAL  WAVE  ATTENUATIONS  IN  THE 
SUBSTRATE  ARE  ASSUMED  TO  REMAIN  CONSTANT  wtTh  DEPTH.  COMPRESSlONAL 
WAVE  ATTENUATION  IN  THE  SEDIMENT  LAYERS  can  vary  with  DEPTH*  AND  CAN 
BE  DISCONTINUOUS  AT  THE  LAYER  INTERFACES. 


MATHEMATICAL  DESCRIPTION 


plmode  REQUIRES  INPUT  FROM  the  node  model  in  order  to  obtain  the 
EIGENVALUES  (KN>»  the  EIGENFUNCTIONS  <UN<Z)),AND  the  attenuation 
terms  in  each  layer  <GNL). 

THE  ATTENUATION  TERMS  ARE  COMPUTED  BY  AN  INTEGRAL  DERIVED  FROM 
FIRST  ORDER  PERTURBATION  THEORY  (SEE  ARL-Tr-bO-13) *  ASSUMING  THAT  THE 
ATTENUATION  IS  CONSTANT  in  EACH  layer.  in  PlMODEi  IT  is  ASSUMED 
that  this  is  true  in  the  water  caw.  the  attenuation  value  in  the 
WATER)  AND  The  SUBSTRATE  <aC  and  AS'  THE  COMpRESSIONAL  AND  SHEAR  WAVE 
ATTENUATIONS*  RESPECTIVELY).  HOWEVER'  ThE  ATTENUATIONS  MAY  VARY 

WITH  depth  In  the  sediment  LAYERS  (AL(Z>>. 

IF  THE  SEDIMENT  ATTENUATIONS  are  constant  in  each  layer  the  mode 
attenuation  for  mooe  number  n  is  COMPUTED  as 

SUM  (  A  *  G  )  ♦  A  •  G  ♦  A*G  ♦  A  *  G 

L  L  Nl  w  NW  C  NC  s  nS 

where  l  Refers  to  a  sediment  layer*  w  the  water*  c  the  compressional 

WAVE  IN  the  SUBSTRATE,  S  ThE  SHEAR  WAVE  IN  THE  SUBSTRATE,  AND  N  THE 
MODE  NUMBER. 

IF  the  sediment  ATTENUATIONS  vary  WITH  Depth,  then  the  mode 
ATTENUATION  FOR  MODE  NUMBER  N,  atn*  is  COMPUTED  as 

SUM  <  D  *  2  *  PI  *  F  /  K  •  INt  >  ♦  A  *  G  ♦  A  *  G  ♦  A  •  G 
L  L  N  L  W  NW  C  Nc  S  NS 

2 

INT  ■  INTEGRAL (  u  <Z)  *  A  (2>  /  C<Z)  > 

L  DZ  N  L 

WHERE  DL  IS  THE  OENSITY  OF  SEDIMENT  LAYER  L*  F  IS  FREQUENCY*  C(Z>  IS 
THE  SOUND  SPEED  AS  FUNCTION  OF  DEPTH,  AND  INTEGRAL (  >DZ  IS  THE 

integral  with  respect  to  depth. 

since  ThE  ATTENUATIONS  Aw,  AC,  AS  AnD  AL(Z)  are  input  in  db/m/kh z 
the  attenuation  sums,  atn*  have  to  be  converted  to  nepers/m 
BEFORE  use  In  the  fielo  CALCULATIONS* 


H  a  AT  •  (F  /  1600.)  *  (LOG  (10.)  /  20.) 

N  N  E 

WHERE  HN  IS  THE  ATTENUATION  OF  THE  NTH  HODF*  ANO  F  IS  THE  FREQUENCY 
(IN  HERTZ). 

THE  FOLLOWING  EQUATION  IS  USED  TO  COMPUTE  THE  COMPLEX  VELOCITY 
POTENTIAL  (VP)  FOR  A  GIVEN  SOURCE  DEPTH  (ZS),  GIVEN  RECEIVER  DEPTH 
(ZR)  AT  A  RANGE#  R» 

I  *  PI  *  Q  *  0  *  SUN  (  H  #  U  (Z  )  •  U  <Z  *  *  EXP  (  -R  *  M  )  ) 

SNNNSNR  n 


(  SL  /  20.  ) 

Q  *  10# 


1 

H  a  H  (  K  *  R) 
NON 


WHERE  DS  IS  THE  DENSITY  OF  THE  LAYER  IN  WHICH  THE  SOURCE  LIES.  SL  IS 
The  SOURCE  LEVEL  <DB  RE  \  MICROPASCAL  at  1  M).  I  IS  THE  COMPLEX  NUMBER 
It  AND  H01(  )  IS  The  HANKEL  FUNCTION  OF  ORnER  ZERO.  FIRST  KINO.  THE 
HANKEL  FUNCTION  IS  APPROXIMATED  BY  THE  FIRST  TERM  OF  THE  ASYMPTOTIC 
expansion, 

I 

H  (K  *  R)  a  SORT (2  ✓  (PI  *  K  •  R) >  *  EXP(  -I  *  (K  *  R  -  PI  /  ♦) >  • 

ON  N  N 

THE  PROPAGATION  LOSS  (PL)  IS  COMPUT'D  BY 


2 

PL  ■  -10.  *  LOG  (  ABS(VP)  )  ♦  SL  . 

10 

The  incoherent  propagation  loss  <ipl>  is  computed  from  the  INCOHERENT 
summation  of  the  mode  terms# 

2 

I  •  PI  •  D  *  SUM  (  (ABS(H  )  *  U  (Z  )  *  u  (z  )  *  EXP  (  -R  #  M  ))  ) 

SN  N  NS  N  R  N 

IPL  •  -lo.  •  LOG  (  IMS  ) 

10 

WHERE  IMS  IS  THE  INCOHERENT  NODE  SUM. 
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the  phase  is  computed  from  the  COMPLEX  Velocity  POTENTIAL* 

-1 

PHASE  *  TAN  (  IM(VP)  /  Rl<VP)  ) 

WHERE  IM(VP)  AND  Rl<VP)  are  the  imaginary  and  REAL  parts* 
respectively*  of  the  velocity  potential. 

the  MAXIMUM  PROPAGATION  LOSS  IS  SET  AT  180  DB.  IF  THE  COMPLEX  FIELD 
IS  ZERO*  THE  PHASE  IS  SET  AT  10**99. 


ON  INPUT 


there  are  two  input  files  required  by  PLMOnE  -  a  formatted  data 
file  supplied  by  THE  USER,  and  The  DISK/TAPE  file  created  by  the 
NORMAL  MOOE  MODEL  NEMESIS. 


lfdisk  (Taped 


This  is  the  disk/tape  file  created  by  the  normal  mode  model*  nemesis 
it  contains  all  the  Relevant  input  data  supplied  to  nemesis  togethek 

WITH  THE  EIGENVALUES,  GROUP  VELOCITIES,  ATTENUATION  TERMS.  AND 
NORMAL  MODES.  SEE  TH£  DESCRIPTION  OF  THE  FORMAT  OF  THIS  FILE  IN 
The  NEMESIS  PROGRAM*  OR  IN  SUBROUTINE  INPMOD  OF  PLMODE. 
this  file  is  read  by  routines  inpmod  and  gftmod. 


LFInP  (TAPES,  INPUT) 


This  is  The  USER-SUPPLIEO  DATA  set  CONTAINING  the  description  of 
The  ATTENUATION  IN  the  ENVlRONMtNT,  and  THF  SOURCE/RECEIVER/RANGE 

geometries. 

this  file  is  read  by  routines  inattn  and  inpsrr. 

the  following  restrictions  are  IMPOSED  On  the  input  TO  PLMODEl 

d)  ALL  DEPTHS  are  REFERENCED  To  The  AlR-WATER  INTERFACE,  AND  NONt 
CAN  fall  IN  The  SUBSTRATE. 

<2)  ATTENUATIONS  MUST  BE  DEFINED  for  EAC"  Layer  that  WAS  SPECIFIED 

in  the  normal  mode  model. 

(3)  FOR  ATTENUATION  PROFILES,  ALL  DEPTHS  MUST  BE  IN  ORDFR  OF 

increasing  depth,  there  must  be  at  least  a  depth-attenuation 
pair  specified  on  the  upper  and  lower  interfaces  of  a  sediment 
layer,  the  last  depth  at  the  bottom  of  a  layer  must  equal  the 
first  depth  at  the  top  cf  ThE  next  layer, 
the  maximum  number  of  depth-attenuation  pairs  is  48. 

(4)  ALL  ATTENUATIONS  MUST  BE  INPUT  IN  DB/meTER/KILOHERTZ. 

<5>  ALL  OEPTHS  MUST  BE  INPUT  IN  METERS*  AMD  MUST  BE  GREATER  THAN  0. 

(6)  ALL  RANGES  MUST  BE  INPUT  IN  KILOMETERS,  AND  MUST  BE  GREATER 
THAN  a. 

(7)  THERE  are  no  RESTRICTIONS  on  The  NUmhfR  OF  SOURCES*  RECEIVERS* 
or  Ranges  that  can  be  input  in  a  data  set  (but  see  -currfnt 
LIMITS-  BELOW) 


INPUT  TO  PlMOOe  IS  VIA  FORMATTED  READ  sTATfM£NTS.  EACH  ITEM  MUST  Bt 
POSITIONED  IN  THE  CORRECT  FIELD  AS  GIVEN  By  THE  FORMAT  STATEMENT. 

CARO  COLUMN  NO 


NO 

OF  CARO 

DESCRIPTION  of  data 

FORMAT 

1 

1-80 

header  information  used  to  describe  data 

baio 

2 

1-10 

attenuation  in  the  water 

Flo. 4 

11-20 

compressional  wave  attenuation  in  the 

SUBSTRATE 

FlO  .4 

21-30 

SHEAR  WAVE  ATTENUATION  In  the  SUBSTRATE 

Flo.4 

IF  SEDIMENT  LAYERS  *ERE  DEFINED  IN  NEMESIS:  CARDS  3,4  ARE  INPUT, 

3  1-A  DESCRIPTION  of  The  TYPE  of  ATTENUATION  Input:  1A4 

COnS  TO  DENOTE  CONSTANT  ATTENUATION  IN  EACH 

sediment  LAYER# 

PROF  TO  DENOTE  DEPTh-aTTENUATION  DESCRIPTION 

for  EACH  SEDIMENT  LAYER 

4A«B, *  FOR  CONSTANT  ATTENUATION  (TONS)* 

attenuation  In  each  sediment  layer 

1-10  ATTENUATION  FOR  FIRST  SEDIMENT  layer 

11-20  ATTENUATION  FOR  SECOND  SEOTMENT  LAYER 

••••  .  •  .  . 

FOR  ATTENUATION  PROFILES  (PRoF>» 

DEPTH-ATTENUATION  PAIRS.  there  must  be  a  pair 
defined  at  the  top  and  bottom  of  each  sediment 


layer 

1-10  DEPTH  Flo. 2 
11-20  SEDIMENT  ATTENUATION  F10.4 
21-30  DEPTH  Flo. 2 
31-40  SEDIMENT  ATTENUATION  FlO. 4 


•  •«.  .  . .  . 

5  BLANK  CARD 

an  ENd-OF-fILE  FOLLOWING  Card  5  INDICATES  the 
end  of  the  data  set  (no  sohrce/Receiver/range 

input  to  COME).  CARDS  1-5-E0F  CAN  BE  REPEATED 
FOR  DIFFERENT  SETS  OF  ATTENTIONS  TO  COMPUTE 
MOOE  ATTENUATIONS  ONLY. 


FlO. 4 
FlO. 4 
FlO. 4 
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l-s  designator  of  the  type  of  output  desired*  ias 


DEPTH 

range 

BOTH 

TO  denote  field,  propagation  loss 

•VS.  DEPTH  FOR  FACH  SOURCE/RANGE 

combination 

TO  DENOTE  field,  PROPAGATION  loss 
.VS.  RANGE  FOR  FACH  SOURCE/RECEIVER 
COMBINATION 

TO  DENOTE  BOTH  TYPES  OF  OUTPUT 

7A,B. • 

SOURCE 

depths 

1-10 

source 

depth 

Flu. 2 

ll-?0 

•  •  • 

source 

•  •  • 

DEPTH 

Flu. 2 

A 

BLANK 

card 

9A«B.  .  RECEIVER  DEPTHS*  8Flw.2 

CAN  BE  INPUT  AS  INDIVIDUAL  DEPTHS  (E.G.* 

10.  1025.),  Op, 

can  be  Input  as  triplets  or  depth  increment 
(preceded  by  a  negative  sirn>*  start  depth  and 

STOP  DEPTH  (E.G*»  -10.  loo.  ?000.>*  OR* 

can  be  input  as  combinations  of  these  two  forms 

OF  INPUT  (E.G.,  1025.  -10.  loO.  2000.  10.  ) 

10  BLANK  CARD 


11  RANGES.  8F10.2 

can  be  Input  in  the  forma i  described 
ABOVE  FOR  the  RLCEIVERS. 

12  BLANK  CARO 

CARDS  6-1?  CAM  BE  repeated  t  INPUT  different 
SOurCE/RECEIvER/RANGE  data  .i  CORRESPONDING 
TO  THt  ATTENUATION  DATA  SET  INPUT  ON  CARDS  1-5. 

an  end-of-fIle  following  card  12  indicates  the 
END  OF  The  entire  data  set. 

THE  NEXT  card  TO  BE  READ  FOLLOWING  AN  ENO-OF-FILE 
HERE  IS  CARD  1. 
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ON  OUTPUT 


LFOUT  (TAPE6) 


this  is  the  printed  output  filf*  it  contains  the  geoacoustic  profile 
describing  the  environment,  the  mode  attenuation  table*  the  complex 
velocity  potential*  COHERENT  and  INCOHERENT  propagation  loss*  and 
PHASE  FOR  EACH  SOURCE/RECEIVER/RANGE  COMBINATION. 

This  FILE  is  WRITTEN  ON  by  ROUTINES  P|_OUT*  PlOUTI*  Pl0uT2»  PLHEAD* 
PLPRNT*  ALERT,  ano  ERROR. 

the  CHARACTER  STRIPS  BOD  AND  EOD  (BEGINNING  AND  END  OF  DATA) 
DELINEATE  THE  START  AND  END  OF  EACH  TABLE  OF  OUTPUT.  nOTE  THAT  WHtN 
THE  PHASE  OF  THE  ACOUSTIC  FIELD  IS  INDETERMINATE  (COMPLEX  FIELD 

equals  zero)*  its  value  is  set  at  io##99.  this  number  exceeds  the 

OUTPUT  field  PROVIDED  FOR  IN  THL  FORMAT*  AND  ON  THE  CYBER  171* 
ASTERISKS  are  OUTPUT. 

The  following  error  messages  may  be  printed  0n  detection  of  an  inpui 

OR  EXECUTION  PROBLEM* 


array  size  diagnostics 

INSUFFICIENT  SPACE  ALLOTTED  To  array  ITITLF  in  main  PROGRAM 

the  value  of  maxtle  in  main  program  Plmode  is  too  Small:  it 

must  BE  GREATER  THAN  OR  EQUAL  THE  NUMBER  OF  WORDS  IN  THE 
FIRST  RECORD  on  the  NEMESIS  output  file. 

insufficient  space  allotted  to  array  buf 

THE  DIMENSION  OF  ARRAY  rUF  IN  SUbROUTinf  INPMOD  (NdUF)  IS 
TOO  SMALL*  it  MUST  be  GREATER  tHaN  OR  EoUaL  TO  The  MAXIMUM 

number  of  layers  specified  in  nemesis  and  the  number  of  words 
in  the  second  Record  on  the  nemesis  disk  file. 

MAX  NO.  UF  LAYERS  ON  MODE  DISK  »GT.  NO.  SET  FOR  THIS  PROGRAM 
THE  VALUE  OF  MAXLAY  IN  MAIN  PROGRAM  PlMuDE  Is  t<jO  SMALL*  IT 
MUST  be  GREATER  Than  OR  EQUAL  to  The  maximum  NUMBER  of  LAYERS 
SPECIFIED  in  nemesis, 

insufficient  space  in  work  array  for  MODES  and  attenuations 

INCREASE  THE  DIMENSION  OF  wOKK  ARRAY  (MaXwRk)  IN  MAIN  PROGRAM 
PLMODE.  SEE  CURRENT  limits  (BELOW)  AS  A  GUIDE  TO  MAXWRK. 

MAX  NO.  OF  SVP  PTS,  On  MODE  DISK  .61.  NO,  SET  FOR  THIS  PROGRAM 
THE  VALUE  OF  MAXSVP  In  MAIN  PROGRAM  Pi_MnOE  IS  TOO  SMALL*  IT 
MUST  BE  GREATER  THAN  OR  EQUAL  TO  THE  MAXIMUM  NUMBER  OF  SOUND 
SPEED  PAIRS  SPECIFIED  In  NFMlSIS. 
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insufficient  work  space  when  allocating  space  for  atten  calcs. 

THE  VALUE  OF  MAXWRk  IN  MAIN  PROGRAM  PlMoDE  IS  TOO  SMALLS 
SEE  CURRENT  LIMITS  (BELOW)  AS  A  GUIDE  To  mAXWRK. 

Insufficient  space  in  attenuation  arrays:  must  be  .ge.  layers-2 

INCREASE  The  DIMENSION  OF  aTTSED.  ATTDEO  ARRAYS  (MAXATT)  in  main 
PROGRAM  PLMODE  TO  AT  LEAST  NUMBER  OF  LAVERS  -  2 

OVERFLOW  in  ATTEN  ARRAY  on  ASSIGNING  ATTen  To  interfaces 

INCREASE  SIZE  OF  ATTSED  AND  ATTDEP  ARRAYS  (MAXATT)  IN  MAIN 
PROGRAM  PLMODE  TO  AT  LEaST  ?*NUM0ER  Of  SEDIMENT  LAYERS 

insufficient  space  in  work  array  for  sedimfnt  atten  calcs 
INCREASE  the  DIMENSION  of  THt  work  ARRAY  (MAXwRK)  In  main 
program  PLMODE.  SEE  CURRENT  LIMITS  (BELOW)  AS  A  GUIDE  TO 
MAXWRK. 


Input  diagnostics 

ERROR  DETECTED  while  reading  from  MODE  FILF 

PROBLEM  WITH  THE  MODE  DISK  FILE:  CHECK  THaT  ThE  FILE  HAS  The 
CORRECT  STRUCTURE*  NUMBER. AND  SIZE  OF  RECORDS. 

EOF  ENCOUNTERED  WHILE  READING  WATER  AND  SUBSTRATE  ATTENUATIONS 

NEGATIVE  WATER  OR  SUBSTRATE  ATTENUATIONS  DETECTED  ON  INPUT 

EOF  ENCOUNTERED  WHILE  READING  SLUlMENT  ATlENUATlON  TYPE 

UNKNOWN  SEDIMENT  ATTENUATION  TYPE  DETECTED  (cOnS  OR  PRoF  ONLY) 

EOF  ENCOUNTERED  WHILE  READING  SLDlMENT  ATTENUATIONS 

negative  sediment  attenuations  detected 

EOE  ENCOUNTERED  WHILE  READING  ATTENUATION  PROFILE 

ZERO  DEPTH  FOUND  IN  ATTENUATION  PROFILE 

IF  BOTH  DEPTH  AND  ATTEN  AHF  ZERO.  THIS  MESSAGE  WILL  NOT  OCCUR !  IT 
IS  EXECUTED  ONLY  IE  DEPTH  IS  ZERO  AND  ATTf N  IS  NON-ZERO. 

OVERFLOW  IN  ATTEN  ARRAYS!  oR»  BLANK  CARD  MUST  FOLLOW  ATTENS 

INCREASE  THE  SIZE  OF  ATTSCn.  ATTDEP  ARRaYS  (MAXATT)  IN  MAIN 
PROGRAM.  OR  REDUCE  The  NUMBER  of  PROFtLF  POINTS  input. 

CHECK  THAT  A  BLANK  CARD  FOLLOWS  THE  ATTENUATION  INPUT 


NEGATIVE  VALUES  DETECTED  IN  ATTENUATION  RpOFILE 
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I.  INTRODUCTION 


V' 

The  normal  mode  model,  NEMESIS,  was  designed  and  implemented  at  the 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin  (ARL:UT), 
to  aid  in  the  investigation  of  low  frequency,  range  invariant,  acoustic 
propagation  effects.  The  mathematical  basis  and  numerical  scheme  of  this 
model  are  described  by  Gonzalez  and  Hawker  (1980) .  NEMESIS  computes  the 
eigenvalues  and  normal  modes  of  the  depth  separated  acoustic  wave  equation. 
Group  velocities  and  modal  attenuation  terms  are  also  computed.  Propaga¬ 
tion  loss  and  velocity  potential  for  a  source  and  receiver  at  a  given 
range  can  be  calculated  from  these  quantities  (Gonzalez  and  Hawker,  1980). 
A  computer  model,  PLMODE,  was  designed  to  perform  these  computations,  and 
the  computer  code  for  both  NEMESIS  and  PLMODE  was  prepared  for  export  to 
other  research  laboratories.  This  report  is  designed  to  accompany  this 
software  and  to  assist  persons  who  wish  to  use  it.  The  contents  include 
descriptions  of  the  modeled  environment,  the  input  data  and  their  format, 
and  the  output  produced  by  the  models.  An  explanation  of  the  diagnostic 
messages  that  may  arise  during  execution  of  the  models  is  also  included. 
Estimates  of  memory  requirements  and  execution  times  are  given  along 
with  sample  data  sets  and  their  resulting  output. 

N 

Although  it  provides  information  beyond  the  scope  of  a  user’s  manual, 
an  appendix,  "Installation  Information,"  has  been  added  to  this  document 
to  aid  in  the  retrieval  of  the  software  from  the  magnetic  tape  on  which 
it  was  written  for  export. 
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MATHEMATICAL  and  physical  description 


NEMESIS  IS  LOW  FREQUENCY*  RANGE" INVARIANT*  Wave  theory  normal  mode 
model*  IT  Was  DESIGNED  to  compute  EIGENVAl  UeS  ano  normal  modes  in  a 
horizontally  stratified  DEEP  ocean  with  single  channel  profiles  ano 
multiple  fluid  sediment  layers  OVERLYING  A  substrate,  sound  speeds 
in  the  WATER  ano  sediment  vary  with  depth*  the  density  is  constant 
within  each  layer*  the  last  layer  is  a  homogeneous*  semi-infinite* 
fluid  or  solid  SUBSTRATE*  l*E.»  the  compressional  and  shear  speeds 

AND  DENSITY  REMAIN  CONSTANT  WITH  DEPTH.  SOUND  SPEEDS  AND  DENSITIES 
Can  BE  DISCONTINUOUS  at  LAYER  INTERFACES.  THE  MODEL  NUMERICALLY 
COMPUTES  THE  OISCRETE  EIGENVALUES  AND  EIGEmFUNcTIONS*  or  normal  modes* 
OF  THE  DEPTH  SePARaTEO  ACOUSTIC  WAVE  EQUATION  rOR  THIS  TYPE  OF 
environment*  THE  CONTINUOUS  PART  OF  the  SPECTRUM  OF  The  solution 
IS  IGNORED*  GROUP  VELOCITIES  ARE  ALSO  COMPUTED* 


DEPTH  EQUATION 


2  2  2  2 
(D  /  DZ  )  u  ♦  (K  <Z)  -  K  )  U  *  0 
N  N  N 


WHERE  U  IS  The  EIGENFUNCTION 
N 

K  is  THE  EIGENVALUE 
N 

K (Z)  *  2*PI*F  /  C<Z> 

C(Z)  IS  THE  SOUND  SPEED 

F  IS  the  source  frequency 


,  ! 


ANY  QUESTIONS  CONCERNING  ThE  APPLICABILITY  OF  THE  MODEL  TO  A 

particular  situation  shoulo  be  directed  to 


KENNETH  E*  HAWKER 

APPLIED  RESEARCH  LABORATORIES 

THE  UNIVERSITY  OF  TEXAS  AT  AUSTIN 

P.  O.  BOX  8029 

AUSTIN*  TEXAS  78712 

TEL*  (512)  835-3407  OR  835-3200 
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ON  INPUT 


the  following  restrictions  are  imposed  on  the  input  to  nemesis 

(1)  there  must  be  at  least  ?  layers  (water  and  substrate)  but  no 
more  THAN  10« 

(2)  ALL  DEPTHS  MUST  BE  IN  INCREASING  ORnER. 

(3)  oepths  at  layer  inteRPacls  must  be  equal, I«e.»  the  last  depth 

at  THE  BOTTOM  OF  A  LAYER  MUST  EQUAL  ThE  FIRST  DEPTH  AT  TOP  OF 
THE  NEXT  LAVER. 

(4)  FIRST  SPECIFIED  DEPTH  MUST  BE  0  (AlR-SURFACE  INTERFACE). 

<5)  THE  SOUND  SPEED  MINIMUM  IN  THE  WATER  AND  SEDIMENTS  MUST  BE  LESS 
than  the  shear  and  cOMPRLSSIONAL  WAVE  SPEEDS  IN  the  substrate. 

(6)  SHEAR  SPEED  OF  ZERO  IN  THE  SUBSTRaTf  IMPLIES  a  FLUID  SUBSTRATE 
AND  A  NONZERO  SHEaR  SPEED  IMPLIfS  a  SOLID  SUbSTraTE.  NO  SHEAR 
SPEEDS  ARt  ALLOWED  IN  ThL  WATER  Or  SEDIMENTS. 

(7)  the  shear  speed  must  be  less  than  tme  compressional  wave  spEto 

IN  THE  SUBSTRATE. 

(8)  IF  A  PRESSURE  RELEASE  hOTTOM  IS  SPECIFIED*  ThfR£  ARE  AN 
INFINITE  number  of  MODES.  The  COMPHrSSlONAL  speed  IS 

USED  AS  THE  PHASE  VELOCITY  CUTOFF.  IN  THIS  CASE,  AND  IS  NOT 
USED  IN  THE  BOUNDARY  CONDITIONS. 

(9)  FOR  SOME  OATA  INPUTS*  DEFAULTS  ARE  RIVEN  IN  PARENTHESES. 

these  can  be  invoked  by  leaving  The  appropriate  data  field 
blank  or  SETTING  it  TO  *. 

(10)  MAXIMUM  NUMBER  OF  DEPTh-SOUND  SPEED  PAIRS  IS  125. 

(ID  MAXIMUM  MODE  *HICH  CAN  b£  COMPUTED  TS  MODE  NUMBER  300. 

(12)  MAXIMUM  NUMBER  OF  M£SH  POINTS  PER  MODE  IS  2000.  THE  NUMBER  OF 
MESH  POINTS  IN  A  LAYER  IS  EQUAL  To  THe  NUMBER  OF  MESH  INTERVALS 
IN  THAT  LAYER  PLUS  1. 

(13)  THE  NUMBER  OF  MESH  INTERVALS  IN  EaCm  LAYER  MUST  BE  EVEN  AND 

greater  than  or  equal  to  there  are  no  mesh  intervals  in  (he 
substrate,  the  solution  is  analytically  computed  in  that 

REGION. 

(14)  ALL  QUANTITIES  must  BE  SPECIFIED  IN  The  mks  system  EXCEPT 
density,  which  IS  Im  ThF  CGS  system. 

(15)  The  OEPTHS  ASSOCIATED  with  both  The  COMPRESSIONAL  and  SHEAR 
speeos  in  the  substrate  are  the  Substrate  intfrface  depths. 

The  COMPRESSIONAL  SPEED  IS  FIRST  &Nr>  THE  SHEAR  SPEED  IS  SECOND. 


INPUT  TO  NEMESIS  IS  VIA  FORMATTED  READ  STATEMENTS.  EACH  ITEM  MUST 
BE  POSITIONED  IN  THE  CORRECT  FIELD  AS  GIVEN  bY  THE  FORMAT  STATEMENT. 

integer  values  and  exponents  must  be  right-justified  in  their  field. 


CARD  COLUMN  no 

NO  OF  CARD  DESCRIPTION  OF  DATA  FORMAT 


1 

2 


3 


4 


1-80 

header  information  used  to  describe  data 

BA10 

1-10 

source  frequency  in  HZ 

Flo. 3 

n-20 

depth  where  the  integrations  are  to  be  hatched* 
must  BE  in  WATER  (DEFAULT  IS  AT  SOUND  SPEED 
minimum  in  thf  water  or  middle  of  the  water 

COLUMN.  DEFAULT  is  RECOMMFNDED) 

Flu.3 

21-30 

wave  NUMBER  k 1 

FlU. 3 

31-40 

wave  NUMBER  K2 

F1U.  3 

41-50 

wave  NUMBER  INCREMENT.  Kl.  K2»  and  INCREMENT 

ARE  FOR  ADDITIONAL  DIAGNOSTIC  INFORMATION.  TH£ 

differential  equation  is  integrated  using  the 
wave  NUMBERS  Kl  TO  K2  IN  The  INCREMENT  SPECIFIED 

FlU. 3 

IN  COLS  41-50.  THE  DISCONTINUITY  IN  THE  DEPTH 

derivatives  of  the  trial  modes  is  tabulated  for 
each  wave  number,  not  Recommended  unless  tracing 
convergence  of  numerical  process,  all  three 
quantities  must  be  posiTivf  to  activate 
DIAGNOSTICS 

1-5  LOWEST  MODE  NUMBER  TO  BE  COMPUTED  (DEFaUlT  IS  1)  15 

6-10  HIGHEST  mode  number  to  Be  COMPUTED  (DEFAULT  IS  15 

TO  COMPUTE  ALL  MODES  UP  TO  MODE  NUMBER  300 > 

OPTIONS  TO  BE  USED  (1*YES*  OsNO> 

5  CORRECT  SOUND  SPEEDS  FoR  CURVATURE  oF  EARTH  15 

10  USE  PRESSURE  RELEASE  SUBSTRATE  15 

15  SUPPRESS  CALCULATION  OF  HIDDEN  DEPTHS*  NOT  15 

recommended 

20  TURN  ON  DIAGNOSTIC  PRINTS  as  MOOES  are  COMPUTED.  15 

not  RECOMMfNDED  UNLESS  TRArlNG  CONVERGENCE  OF 
NUMERICAL  PROCESS 


5 


1-10 

11-15 


DENSITY  OF  LAYER 

NUMBER  OF  MESH  INTERVALS  In  LAYER  (DEFAULT  IS 
maximum  ALLOWABLE) 


Flu. 3 
IS 


6A«B« • 

SOUND 

SPEED  PROFILE  IN  THIS  LAYER 

1-10 

OEPTH 

11-20 

SOUND 

SPEED 

21-30 

DEPTH 

•  •  f 

•  •  • 

71-80 

SOUND 

SPEED 

7 

BLANK 

CARD 

ft 

9 


1-10 

1-10 

11-20 

21-30 

31-40 


REPEAT  CARDS  5-7  FOR  EACH  SEDIMENT  LAYER 

density  of  substrate 

DEPTH  OF  SEOImENT/SUBSTRaTf-  INTERFACE 
COmPRESSIOnAL  speed 

depth  of  SEDImENT/SUBSTRaTf  INTERFACE 

shear  speed 


OK  OUTPUT 


NEMESIS  USES  FILES  LFOUT,  lTAPE1»  AND  LFDIaG  as  output  files 
ON  LFOUT 


NEMESIS  USES  FILE  LFOUT  to  OUTPUT  ALL  PRINTED  INFORMATION  SUCH  AS 
INPUTS  AND  QUANTITIES  TNaT  WERE  CALCULATED  In  THE  PROGRAM  SUCH  AS 
EIGENVALUES*  GROUP  VELOCITIES,  AND  ATTENUATION  TERMS.  the 
ATTENUATION  TERMS  FOR  EACH  MODE  are  LISTED  IN  the  FOLLOWING  ORDER: 

ATTENUATION  in  water  LAVER 
ATTENUATION  in  first  sediment  layer 

ATTENUATION  in  LAST  SEDIMENT  LAYER 
ATTENUATION  in  SUBSTRATE  (COMPRESSIOmAl) 

ATTENUATION  in  SUBSTRATE  (SHEAR) 

ON  LED I AG 


NEMESIS  USES  LFDIAG  FOR  PRINTED  DIAGNOSTIC  MfSSAGES.  CAUTION  -  IF 
DIAGNOSTIC  PRINT  OPTION  IS  TURNED  ON#  THIS  file  can  be  VERY  LARGE. 

The  FOLLOWING  ERROR  ANO  OIaGNOSHC  MESSAGES  aRE  generated  by  NEMESIS: 

FATAL  INPUT  ERRORS 

only  io  layers  are  allowed 

at  least  two  layers  MUST  be  USED  TO  Run  NEMESIS 
SOUNO  SPEED  ARRAY  OVERFLOW 
all  depths  MUST  BE  IN  INCREASING  order 
depths  at  layer  INTERFACES  must  be  EQUAt 

each  layer  must  have  AT  LEAST  TWO  DEPTH-SOUND  SPEFD  points 
depth  AT  AIR-WATER  INTERFACE  must  be  zero 


SOUND  SPEED  MINIMUM  OF  THE  WATER  AND  SEDIMENT  LAYERS 

must  he  less  than  the  compressional  velocity 

SHEAR  SPEED  MUST  Be  LESS  Than  The  CoMPWfSsIONAL  speed 

if  THE  SHEAR  SPeEO  IS  NONZFRO,  THE  SOtJNn  SPEED  MINIMUM 
must  Be  LESS  THAN  THE  SHEAR  SPEED 

OPTIONS  MUST  BE  SPECIFIED  wlTh  A  0  OR  1 

TOO  MANY  SPEEDS  FOR  SURSThaTl 

fatal  errors  (numerical  and  programming) 

INTEGRATION  COUNT  exceeds  maximum  ALLOWABLE*  mooe  no  S 

all  MATCH  POINTS  EXHAUSTED,  CONVERGENCE  NoT  ACHIEVED,  FOR  MODE 

ERROR  IN  NUMROV,  CALLED  FROM  INTGR1 *  IER  = 

error  in  subroutine  INTRP  after  CALL  TO  IcSlCU 

ERROR  IN  SUBROUTINE  INTrP  after  CALL  to  IcSEVU 

ERROR  in  CALL  To  Secant  FOR  MODE  NUmBEW. 

ERROR  IN  COMPUTATION  of  RAYLlIGH  VELOCITY 

ERRORS  DETECTED  IN  IMSL  ROUTiNLS 

***  TERMINAL  ERROR,  IER  *  FROM  IMSl  ROUTINE 

***  WARNING  WITH  Fix  ERROR  IER  =  From  jMSL  ROUTINE 

***  WARNING  ERROR  IER  a  FROM  IMSL  ROUTINE 

•**  undefined  error  ier  =  from  imSl  routine 


diagnostic  prints  for  tracing  convergence  of  numerical  process 
(for  more  detailed  Information  see  Gonzalez  and  hawker#  i960) 


THE  CURRENT  MODE  NUMBER  BEING  COMPUTED  1 
MODENO  s 

ONE  OF  THE  FOLLOWING  IS  PRINTED  FOR  EACm  TRIAL  EIGENVALUE: 

FIRST  ESTIMATE 
second  estimate 
secant  iterate 

right  BOUNO  ESTIMATE  FROM  refine 
left  BOUND  ESTIMATE  FROM  refine 
right  BOUND  ESTIMATE  from  zeroin 
left  BOUND  ESTIMATE  from  ZEROIN 
BISECTION  ESTIMATE 

THE  TRIAL  EIGENVALUE*  NUMBER  OF  ZEROS  Iki  the  TRIAL  EIGENFUNCTION# 

and  integration  NUMBER  I 

KN  r  NZEROS  ®  INTEGRATION  NO  s 

THE  FIVE  POINTS  CENTERED  AT  THE  MATCH  PolNT  FOR  EACH  UPPER  AND 

lower  half  of  trial  mode* 

U-UPPER  U-LOWER 

the  depth  derivative  EVALUATED  at  ThE  m#TCH  point  For  each 
UPPER  AND  LOWER  HALF  OF  TRIAL  MODE: 

UPPER  DERV  s  LOWER  OERV  = 

ThE  DISCONTINUITY  In  thf  depth  DERIVATIVES  at  the  match  point* 
I.E.#  DISC  a  UPPER  OERV  -  LOWER  DERV.  RaTIO  IS  THE  SCALING 
FACTOR  that  WAS  APPLIED  to  the  UPPER  HAI  F  To  FORCE  CONTINUITY 
IN  THE  UPPER  AND  LOWER  HALVES  OF  TRIAL  mOdEJ 

DISCONTINUITY  *  RATIO  a 


} 

i 

I; 

i 

i 


! 

i 
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DIAGNOSTIC  PRINTS  FOR  TRACKING  DISCONTINUITY  (PhI)  IN  first 
DERIVATIVE  OF  TRIAL  EIGENFUNCTIONS  for  a  range  of  trial 
EIGENVALUES  for  a  FIXED  MATCH  point: 

PHI  ZEROS  ZMATCH 


ON  LTAPE1 


NEMESIS  GENERATES  A  DISK/TAPE  FILE  ON  lTAPFI  ThAT  CONTAINS  ALL  OF  THE 
RELEVANT  INPUT  DATA  TOGETHER  WITH  ALL  OF  T«E  EIGENVALUES,  GROUP  VELOC¬ 
ITIES*  ATTENUATION  TERMS,  aND  NORMAL  MOOES.  NOTE  -  WITH  THE 
EXCEPTION  of  The  HEAOER  INFORMATION,  all  Ut'ANTiTlES  are  OUTPUT  AS 
FLOATING  POINT  NUMBERS.  LTaPEI  has  the  FOLlOwlNG  STRUCTURE* 


LOGICAL  WORD  NO  DESCRIPTION  OF  RECORD  NO  OF  WU*0S 

RECORD  NO  IN  RECORD  PER  RECORD 


1 

1-8 

header  information 

8 

2 

1 

NUMdER  OF  LAYERS 

20 

2 

maximum  allowable  layeps 

3 

0.0 

4 

frequency  in  hertz 

5 

COMPResSlONAL  SPEED  IN  SUBSTRATE 

6 

shear  speed  in  substrate 

7 

Rayleigh  spled  at  substrate  interface 

8 

pressure  release  substrate  (isYES*o*NO) 

9 

last  modl  computed 

10 

FIRST  MODE  COMPUTED 

11-20 

O.o 

3 

1-10 

Index  of  thl  last  mesh  point  in  each  layer 
of  the  modes,  beginning  with  layer  i 

NO 

4 

1-10 

index  of  thc  last  point  in  each  layer  of 

The  input  sound  speed  PROFILE 

NO 

5 

1-10 

depth  at  each  layer  interface 

NO 

6 

1-10 

mesh  InCPpmLNT  (In  MEIFRs)  FOR  EACH  LAYER 

NO 

7 

1-10 

density  of  lach  layer 

NO 

ir-n 


8 

l-Nl 

DEPTHS  CORRESPONDING  To  THE  INPUT  SOUND 
SPEEO  PROFILE 

Nl 

9 

1-N1 

SOUNO  SPEEDS  CORRESPONMNG  TO  THE  INPUT 
SOUND  SPEED  PROFILE 

Nl 

10  - 
(N3«9) 

1-N2 

NORMAL  MODES  EVALUATED  AT  EACH  mesh  POINT 

(there  is  one  record  pfr  mode) 

N2 

N3M0 

1-N3 

EIGENVALUES  CORRESPONDING  TO  EACH  mode 

N3 

N3.J1 

1-N3 

GROUP  VELOCITIES  OF  EACH  MODE 

N3 

N3.12 

1-N5 

attenuation  terms  per  i ayer  per  mode 

attenuation  terms  FOR  FACH  mode  in  each 

NA 

layer*  beginning  with  i  ayer  i.  both  shear 
And  COMPRESSIONAL  ATTENUATIONS  are  computed, 
therefore*  there  are  ac  many  attfnuation 
terms  for  each  mode  as  there  are  layers*i 

(N5.D-2N5  attenuation  terms  for  second  mode  computed 

<2N5^1)-3NS  ATTENUATION  TERMS  FOR  THIRD  MODE  COMPUTED 


( <N3-1)*N5*1>  ATTENUATION  TERMS  FOR  |  AST  MODE  COMPUTED 
-N3*N5 

end  -  of  -  file 


total  NUMBER  of  words  IN  FILL  LTAPEl  IS 
28  ♦  5*N0  ♦  2*N1  ♦  (2*N?*N5) *N3 


DEFINITIONS 

no  MAXIMUM  ALLOWABLE  LAYERS 

N1  TOTAL  NO.  OF  POINTS  IN  INPUT  SOUND  SPEED  PROFILE.  ACCESS 
WORD  NUMBER  (LAYERS-1)  IN  *Th  RECORD 
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total  NO.  OF  mesh  POINTS  in  mooes.  ACCESS  word  number 
(LAYERS-!)  in  3RD  RECORD 

TOTAL  NO.  OF  MOOES.  SUBTRACT  FIRST  MODE  COMPUTED  FROM  LAST 
MODE  COMPUTEO  AND  ADO  1 

N3  •  NUMBER  OF  LAYERS 

number  OF  LAVERS  ♦  1 


required  RESIDENT  routines 


STANDARD  FORTRAN 

ABS»AL0610*AMAXl*AMlNl«FLOAT«lABS»lFlX,MAXO»MINOtMODtSQRT 


NONSTANDARD  -  THESE  HUST  BE  WRITTEN  AT  EaCm  nEW  INSTALLATION 
EOF  (U) 

purpose  function  that  tests  for  an  end  of  file 

CONDITION  FOLLOWING  a  READ 
FORMAL  PARAMETERS  U  -  UNIT  (FILE  name) 
value  RETURNED  ZERO  IF  NO  END  OF  FILE  IS  ENCOUNTERED 

NONZERO  otherwise 

IBUFOUT  (U» IFRhAT* ARRAY*NPTS) 

purpose  function  that  buffers  out  a  contiguous 

BLOCK  OF  WORDS 

FORMAL  PARAMETERS  U  -  UNIT  (FIlE  NAME) 

IFRMAT  -  ONE  for  BINARY  FORMAT 

ARRAY  -  ADDRESS  oF  FIRST  WoRD  IN  BLOCK 

npts  -  number  of  words  to  be  buffered  our 
VALUE  RETURNED  zero  IF  SUCCESSFUL 

NONZERO  OTHERWISE 


current  limits 


VARIABLE 

DESCRIPTION 

LIMIT 

meshmx 

maximum 

number  of  mesh  PoiwTS 

2000 

lyrmax 

MAXIMUM 

number  of  layers 

io 

mxinpt 

MAXIMUM 

NUMBER  of  DEPTH-SOUND 

SPEED  PAIRS 

WHICH  CAN  BE  INPUT 

nmoomx 

MAXIMUM 

mode  number  which  CAN 

be  computed 

300 

array  DIMENSIONS 


ARRAY 

Array  DIMENSION 

AKSQK 

meshmx 

depthi 

MeSHMX 

speed 

MeShMX 

U 

MeShMX 

aksqrh 

lyrmax 

buffi 

LyRmAX 

BUFF2 

LYRMAX 

BUFF3 

lyrmax 

chalf 

LYRmAX 

delz 

lyrmax 

H 

LYRMAX 

MESH 

LYRmAX 

net 

LYRMAX 

NPTS 

LYRmAX 

rho 

LYRMAX 

ATTEN 

I yrmax*i 

z 

mxinpt 

c 

MXINPT 

COfcF 

MXINPT*3 

bufakn 

nmOdmx 

bufgv 

NmODMX 

BUFATT 

NMODMX* (lYRMAX *1 ) 

IHEADR 

BO  CHARACTERS 

COMPUTER  RESOURCES 


NEMESIS  NAS  PROGRAMMED  ON  A  CONTROL  DATA  CORPORATION  CYBER  171* 

under  the  nos  1.3  level  ass  operating  systfm.  this  computer  has 

A  60-BIT  WORD  WITH  10  CHARACTERS  PER  WORD.  THE  FORTRAN  COMPILER 
is  fortran  four  extended  u.7#4B5>  which  is  described  in  cdc 
publication  60497000*  revision  d.  efforts  were  made  to  maintain 

ANSI  STANDARD  FORTRAN  CODE.  PORTIONS  OF  TwE  CODE  ARE  UNAVOIDABLY 

Installation  dependent#  such  as  the  number  of  characters  per  word. 

ON  THIS  MACHINE#  THE  LAST  wORO  address  (LWA)  of  NEMESIS  was 
62121  OCTAL  WOROS.  ThIS  caN  BE  REDUCED  RY  DfCREaSIng  THE  SIZES 
OF  THE  ARRAYS  IN  THE  MAIN  PROGRAM  AND  APPROPRIATE  COMMON  BLOCKS# 

AS  DESCRIBED  IN  THE  SECTION#  CURRENT  LIMITS. 

nemesis  was  run  with  a  typical  oeep  water  profile  with  two 
sediment  layers  and  a  fluio  substrate,  eigenvalues,  eigenfunctions, 

GROUP  VELOCITIES#  AND  ATTENUATION  TERMS  WERE  COMPUTED  AND  STORED  ON  A 
FILE#  THE  TOTAL  EXECUTION  TIME  FOR  THE  14o  mOoES  WITH  1999  MESH 
POINTS  PER  MOOE  WAS  194.15  CP  StCONDS  (OPTIMIZED).  SEE  ARL-TM-BO-6# 
CHAPTER  3. FOR  FURTHER  DETAILS  OF  THIS  EXAMPLE# 


III.  NEMESIS:  EXAMPLES 


Two  data  sets  and  the  corresponding  printed  output  from  NEMESIS  are 
included  in  this  chapter.  The  modeled  environment  is  shown  in  Fig.  1. 

It  consists  of  four  layers:  water,  two  sediment  layers,  and  the  substrate. 

3 

The  density  in  each  of  the  four  layers  is  1.0,  1.5,  1.75,  and  2.0  g/cm  , 
respectively.  The  sound  speeds  in  the  sediment  layers  vary  linearly  with 
depth,  and  there  is  a  discontinuity  in  sound  speed  at  the  water- sediment 
interface.  Note  that  the  sound  speed  must  be  specified  at  the  top  and 
bottom  interface  of-  each  water  and  sediment  layer,  and  that  within  the 
substrate  the  shear  and  compressional  speeds  must  be  constant. 

The  first  data  set  does  not  include  shear  waves  in  the  substrate; 
i.e.,  the  substrate  is  a  fluid,  and  thus  the  shear  speed  is  zero.  At 
50  Hz  there  are  149  modes,  as  shown  in  the  printed  output.  The  five 
attenuation  terms  listed  for  each  mode  pertain  to  the  compressional  wave 
in  the  water,  the  sediment  layers  and  the  substrate,  and  the  shear  wave 
in  the  substrate,  in  that  order.  Because  shear  waves  are  not  included, 
the  shear  wave  attenuation  term  is  identically  zero  for  all  modes.  In 
this  data  set  the  option  to  suppress  hidden  depth  calculations  is  not 
chosen.  Thus,  for  modes  with  turning  points  high  above  the  water-sediment 
interface  (modes  1  to  49),  the  sediment  attenuation  terms  are  zero.  Note 
that  the  phase  velocities  of  the  modes  are  bounded  by  the  minimum  sound 
speed  in  the  water  and  the  substrate  compressional  speed.  This  computa¬ 
tion  requires  194  cp  sec  on  a  Control  Data  Corporation  CYBER  171  when 
NEMESIS  is  compiled  under  maximum  optimization. 

The  second  data  set  describes  the  environment  in  Fig.  1  with  a  solid 
substrate.  The  substrate  compressional  speed  is  5000  m/sec  and  the  shear 
speed  is  2700  m/sec.  At  50  Hz  there  are  282  modes  with  phase  velocities 
bracketed  between  the  minimum  sound  speed  in  the  water  and  the  shear 
speed.  This  computation  requires  347  cp  sec  on  the  CYBER  171  in  the 
optimized  mode. 


1 


stitutes  the  data  set  of  the  first  example 
:1  NEMESIS,  for  this  report. 
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Ml'NK  PRO  mi.  Cl  -  1500,  Z1  ■  1000,  IPS  -  .0074,  S  -  1233.33 
50. 


1.000 

0.000 

1527.081 

51.000 

1523.539 

101.000 

1520.411 

152.000 

1517.543 

202.000 

1515.024 

253.000 

1512.728 

303.000 

1510.725 

354.000 

1 508 .9 15 

404.000 

1507.351 

455.000 

1505.952 

505.000 

1504.760 

556.000 

1503.712 

606.000 

1502.836 

657.000 

1502.065 

707.000 

1501 .477 

758.000 

1500.978 

808.000 

15e8.598 

859.000 

15e0.314 

509.000 

15e0.127 

960.000 

1500.024 

1010.000 

1500.001 

1061.000 

1500.053 

1111.000 

1500.170 

1162.000 

1500.352 

1212.000 

1500.587 

1263.000 

1500.880 

1313.000 

1501 .216 

1364.000 

1501 .6es 

1414.000 

1502.024 

1465.000 

1502.492 

1515.000 

1502 .935 

1566.000 

1503.521 

1616.000 

1504.076 

1667. 0e« 

1504.669 

1717.000 

1505  .276 

1766.000 

1505.919 

1818.000 

1506.570 

1869.000 

1507.254 

1919.000 

1507  .943 

1970.000 

1508.662 

2020.000 

1505.383 

2071.000 

1510.133 

2121.000 

l5ie.S80 

21'’2.000 

1511.655 

2222.000 

1512.426 

2273.000 

1513.223 

2323.000 

1514.013 

2374.000 

1514.828 

2424.000 

1515.635 

2475.000 

1516.465 

2525.000 

1517  .286 

2576.000 

1518.130 

2626.000 

1518.563 

2677.000 

1519.618 

2727.0.0J! 

1520.561 

27*8.000 

1521.525 

2828.000 

1522.377 

2879. 00e 

1523.249 

2929.000 

1524.108 

2980.000 

1524.988 

3030.000 

1525.853 

3081.000 

1526.738 

3131.000 

1527  .606 

3182.000 

1528.499 

3232.000 

1529.374 

3283.000 

1530.266 

3333.000 

1531.147 

3384. e00 

1532.045 

3434.000 

1532.926 

3485.000 

1533.827 

3535.000 

1534  .712 

3586.000 

1535.616 

3636.000 

1536.503 

3667.000 

1537.408 

3737.000 

1538.297 

3786.000 

1539.205 

3638.000 

1540.095 

3669.000 

1541.005 

3539.000 

1541.697 

3996  .000 

1542.807 

4040.000 

1543.700 

4091.000 

1544.612 

4141.000 

1545.508 

4192.000 

1546.419 

4242.000 

1547.314 

4293.000 

1548. 22* 

4343 . 000 

1549.123 

4394.000 

1550.037 

4444.000 

1550.934 

4495. 0C0 

1551  .849 

4545.000 

1552.746 

4596.000 

1553.661 

4646.000 

1554.558 

5000.000 

1560.917 

1.500 

5000.000 

1550. ee0 

5100.000 

1650.000 

1.750 

5100.000 

1650.000 

5200.000 

1680.000 

2.000 

5200.000  1700.000 


B.  PRINTED  OUTPUT  (1) 


The  following  printout  is  the  output  generated  by  data  set  (1)  of 
NEMESIS. 


NEWS  IS  ARL-UT  NORMAL  MODI  MODEL 

VERSION  1.2 


MUNE  PROFILE,  Cl  -  1500,  21  -  1000,  EPS  «  .0074,  B  «  1233.33 


FREQUENCE  (HZ) 


50.000 


NUMBER  OP  EXISTING  MODES 
FIRST  MODE  COMPUTED 
LAST  MODE  COMPUTED 


146 

1 

149 


DEPTH  (M)  AT  WHICH  MODES  ARE  MATCHED 


< 


1010.000 


NO  OPTIONS  SELECTED 


GEOACOUSTIC  PROFILE 


LATER 


1 


TOTAL  LOCAL  COMPRESS IONAL  SHEAR 


DEPTH 

(K) 

DEPTH 

(M) 

0 .000 

0.000 

51.000 

51.000 

101.000 

101.000 

152.000 

152.000 

202.000 

202.000 

253.000 

263.000 

303.000 

3e3.000 

354.000 

354.000 

404 .000 

404.000 

455. cec 

455.000 

505.200 

505.000 

556.000 

556.000 

606.000 

606.000 

657.C00 

657.000 

707.000 

707.000 

758.000 

758.000 

606.000 

808.000 

859.000 

859.000 

909. 0e0 

909.000 

960.000 

960.000 

1010.000 

1010.000 

1061.000 

1061.000 

1111.000 

1111.000 

1162.000 

1162.000 

1212.000 

1212.000 

1263.000 

1263.000 

1313.000 

1313.000 

1364.000 

1364.000 

1414.000 

1414.000 

1465.000 

1465.000 

1515.000 

1515.000 

1566.000 

1566.000 

1616.000 

1616.000 

1667.000 

1667. 0ee 

1717.000 

1717.000 

1768.000 

1768.000 

1816.020 

lsie.eee 

1869.000 

ie69.000 

1919.020 

1919.000 

1970.000 

1970 .tee 

2020.000 

2020.000 

2071  .000 

2071.000 

2121 .000 

2121. eee 

2172.000 

2172.000 

2222.000 

2222.000 

2273.000 

2273.000 

bAVE  SPIES  bAVE  SPEED 

(M/SEC)  (M/SEC) 


1527. Col 
1523.539 
1520.411 
1517.543 
1515.024 
1512.728 
151? .725 
1508.915 

1507.351 
1505.952 
1504.760 
1503.712 
1502.636 
1502.085 
1501 .477 
15e0.978 
1500.598 
1500.314 
1500.127 
1500.024 
1500.001 
1500.053 
1500.170 

1500.352 
1500.587 
1500.860 
1501.216 
1501.603 
1502.024 
1502.492 
1502.985 
1503.521 
1504.076 
1504.669 
1505.276 
1505.919 
1506.570 
1507.254 
1507.943 
1508.662 
1509.383 
1510.133 
1510 .880 
1511.655 
1512.426 
1513.223 


RATLEIOH  SAMPLING 

WAVE  SPEEI  DENSITY  POINTS 

th/SEC)  IG/'CC) 


1.020  1921 


2374.000 

2374.000 

1514.828 

2424.000 

2424.000 

1515.635 

2475.000 

2475.000 

1516.465 

2525.000 

2525.000 

1517.286 

2576.000 

2576.000 

1518.130 

2626.000 

2626.000 

1518.963 

2677. e00 

2677.000 

1519.818 

2727.000 

2727.000 

1520.661 

2778.000 

2778.000 

1521.525 

2828. eee 

2828.000 

1522.377 

2879.000 

2879.000 

1523.249 

292  9.000 

2929.000 

1524.108 

2980.000 

2960. eee 

1524.988 

3030.000 

3030.000 

1525.  B53 

3081.000 

3081.000 

1526.738 

3131.000 

3131.000 

1527.608 

3182.000 

3182.000 

1528.499 

3232.000 

3232.000 

1529.374 

3283.000 

3283.000 

1530.268 

3333.000 

3333.000 

1531 .147 

3364.000 

3384.000 

1532.045 

3434.000 

3434.000 

1532.926 

3465.000 

3485.000 

1533.827 

3535.000 

3535.000 

1534.712 

3586.000 

3566.000 

1535.616 

3636.000 

3636.000 

1536.503 

3687.000 

3687.000 

1537.408 

3737.000 

3737.000 

1538.297 

376e.000 

3788.000 

1539.205 

3838.000 

3638.000 

1540.095 

3889.000 

3989. eee 

1541.005 

3939.000 

3939.000 

1541.697 

3990.000 

3990.000 

1542.807 

4040.000 

4040.000 

1543.700 

4091.000 

4 091.000 

1544.612 

4141.000 

4141.000 

1545.506 

4192.000 

4192.000 

1546.419 

4242.000 

4242.000 

1547.314 

4293.000 

4293.000 

1548.227 

4343.000 

4343.000 

1549.123 

4394.000 

4394 . 000 

1550 .037 

4444 .000 

4444.000 

1550.934 

4495. 0e0 

44S5.000 

1551.849 

4545.000 

4545.000 

1552.746 

4596.000 

4596.000 

1553.661 

4646.000 

4646. 0e0 

1554.558 

5000.000 

5000.000 

1560.917 

4 

5000.000 

0.000 

1550 .000 

1  •5*2 

39 

_  _  _ 

5100.000 

100.000 

1650.000 

3 

5100. e00 

0.000 

1650.000 

1.750 

36 

_  _  _ 

5200.000 

If  '00 

1680.000 

4 

5200.000 

„  .000 

1700.000 

e  .000 

0.000 

2.0<.V 

t 
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MODE  CHARACTERISTICS 


i 

j 


MODE 

WAVE  NUMBER 

PHASE  VELOCITT 

CROUP  VELOCITT 

ATTENUATION  TERMS 

NUMBER 

(1/M) 

(M/SEC) 

(M/SEC) 

1 

.203369840 

1500.4991 

1500.0009 

.1000166031*01 

0. 

0. 

0. 

0. 

2 

.209231409 

1501.4919 

1500.0084 

.1000494291*01 

0. 

0. 

0. 

0. 

3 

.209094015 

1502.4785 

1500.0227 

.ieezoi806E*ei 

0. 

0. 

0. 

0. 

4 

.208957603 

1503.4594 

1500.0443 

.100113733E+01 

0. 

0. 

0. 

0. 

5 

.208822191 

1504.4343 

1500.0725 

.1001452231*01 

0. 

0. 

0. 

0. 

6 

.208687766 

1505.4034 

1500. ie74 

. 1001762 82E*01 

e. 

0. 

0. 

0. 

7 

.208554294 

1506.3668 

1500.1480 

. 1 0020 6944 E* 01 
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.186577176 

1683.6033 

1410.8728 

.940461261 £♦00 
.3974096011-01 
.9987731 52E-01 
.1182 57902 E-01 
0. 

|  144 

1 

1  ' 

! 

.186321963 

1686.1094 

1403.2262 

.9714219861+00 
.  36223 170eE-ei 
.7593564861-01 
.1211497451-01 

e. 

) 


i 

!  • 

? 

I 

i 

» 


•i 


v>  •*  :  -v 


145 

.186056320 

1688.5170 

1395.2417 

.1003569881*01 

.3010335e5E-01 

.5436231551-01 

.1165410781-01 

0. 

146 

.185781564 

1691.0142 

1369.3634 

.1024667361*01 
.2516534531-01 
.411 57e844E-01 
.1166977391-01 
0. 

147 

.185500676 

1693.5748 

1365. 7072 

.1034436501*01 
.2215966551-01 
.3*8244351 1-01 
.1 37^36973 E-Cl 
0. 

148 

.185216616 

1696.1721 

1364.4528 

.1 C3257264E*01 
.20S'»6  54051-01 
.3351949571-01 
.1954800181-01 
0. 

149 

.184933832 

1696.7656 

1369.4109 

.1003922471*01 

.2110512771-01 

.3692474731-01 

.4330241471-01 

0. 


The  following  printout  constitutes  the  data  set  for  the  second 
example  generated  by  NEMESIS. 


MCKK  PROrilE 
50. 


1.000 
e.eee 
202.000 
404.000 
606.000 
606.000 
1010.000 
1212.000 
1414.000 
1616.000 
1615.000 
2020 . 000 
2222.000 
2424.000 
2626 .000 
2829.000 
3030.000 
3232. 00* 
3434.000 
3636.000 
3P35.000 
4040.000 
4242.000 
4444.000 
4646.000 

1.500 
50t  2  .020 

1.750 

5100.000 

2.000 

5200.000 


,  Cl  -  1500,  Z1  -  1000,  EPS  -  .0074,  £  *  1233.33 


1527.081 
1515.024 
1507.351 
1522.836 
1502.596 
1500. 001 
1500.587 
1502.024 
1504.076 
1506.570 
1529.383 
1512.426 
1515.635 
1518.963 
1522.377 
1525.853 
1529.374 
1532.926 
1536.503 
1540.095 
1543.700 
154 ’'.314 
1550.934 
1554. 55P 


51.000 

253.000 

455.000 

657.000 

859.000 

1061.000 

1263.000 

1465.000 

1667.000 

1669.000 

2071.200 
2273.000 
24’75.000 
2677.000 

2879.200 
3061.00.0 
3283.000 
3485.000 
3687.000 
3?86 . 000 
4091.000 
4293.000 
4495.000 
5002.000 


1523.539 
1512.728 
1505.952 
1502.085 
1500.314 
1500.053 
1500.680 
1502.492 
1504.669 
1507.254 
1510.133 
1513.223 
1516.465 
1519.816 
1523.249 
1526 .73e 
1530.268 
1533.827 
1537  .40e 
1541 .005 
1544.612 

1548.227 

1551.849 

1560.917 


101.000 
303.000 
505.000 
707.000 
909.000 
1111.000 
1313.000 
1515.000 
1717.000 
1919.00? 
2121.000 
2323. 040 
2525.000 
2727.0a? 
2929.000 
3131. 020 
3333.000 
3535. 000 
3737. ecfc 
3939.002 
4141.000 

4343. ?ee 

4545.000 


1520.411 
1510.725 
1504.760 
1501.477 
1500.127 
1500.170 
1501.216 
1502.985 
1505.276 
1507.943 
1510.880 
1514.013 
1517  .286 
1520.651 
1524  .103 
1527  .50? 
1531.147 
1534  .712 
1538.297 
1541.897 
1545.506 
1549.123 
1552.746 


152.000 
354.000 
556.000 
758. e00 
960.000 
1162.000 
1364. 0e0 
1566.000 

l?6e.000 

1970.000 
2179.000 
2374.000 
2576.000 
277 t.000 
2980.000 
3152.000 
3364.000 
3586.000 
3768. e00 
3660.000 
4192.000 
4394. Z0Z 
4596. 0Z0 


1517.543 
1508.915 
1503.712 
1500.979 
1500.024 
1500.352 
1501.603 
1503.521 
1505.919 
1505.662 
1511.655 
1514.828 
1513.130 
1521.525 
1524.986 
1528.499 
1532.045 
1535.616 
1536. 2e5 
1542.807 
1546.419 
1550.037 
1553.661 


1550.000  5100.00?  lese.eee 
1650. ?ee  5200.000  ieee.000 
500Z.00Z  5200.000  z’00.000 
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PRINTED  OUTPUT  (2) 


The  following  printout  is  the  output  of  data  set  (2)  from  NEMESIS. 


nimes is  arl-ut  normal  moll  model 

VERSION  1.2 


MUNI  PROFILE,  Cl  «  1500,  Z1  *  1302,  EPS  *  .0074,  B  -  1233.33 


FREQUENCY  (EZ)  50  .e00 

NUMBER  or  EXISTING  MODES  292 

FIRST  MODE  COMPUTED  1 

LAST  MODE  COMPUTED  292 

DEPTH  (M)  AT  VHICfi  MODES  ARE  MATCHED  1010.302 


NO  OPTIONS  SELECTED 
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- -  ' 


GEOACOUSTIC  FROIILE 


TOTAL  LOCAL  COMPRESS IONAL  SHEAR  RATLEIGH 

LATER  LEPTH  PEPTH  WAVE  SPEEP  WAVE  SPEFE  WAVE  SPEEE  PENS  ITT 

(M)  (M)  (M/SEC)  (M/SEC)  (M/SEC)  (G/CC) 


0.020 

0.000 

1527.081 

SI .000 

51.000 

1523.539 

lei  .eez 

101.000 

1520.411 

152.000 

152.000 

1517.543 

202.000 

202.000 

1515.024 

253.000 

253.000 

1512.728 

303.000 

303.000 

1510.725 

354.000 

354.003 

1508.915 

404.000 

404.000 

1507.351 

455.000 

455.000 

1505.952 

505.000 

505.000 

1504 .760 

556.000 

556.000 

1503.712 

606.000 

606. e00 

1502. e36 

657.000 

657.000  • 

1502.085 

707.000 

707.000 

1501.477 

756.000 

756.000 

1500.978 

606.000 

806.000 

1500.598 

659.000 

859.000 

1500.314 

909.000 

909.000 

1500.127 

960.000 

960.000 

1500.024 

1010.000 

1010.000 

1500.001 

1061.000 

1061.000 

1500.053 

1111.000 

1111.000 

1500.170 

1162.000 

1162.000 

1500.352 

1212.000 

1212.009 

1500.587 

1263.000 

1263.000 

1500.880 

1313.000 

1313.000 

1501.216 

1364.000 

1364.000 

1501.603 

1414.000 

1414.000 

1502.024 

1465.000 

1465. oee 

1502.492 

1515.000 

1515.000 

1502.985 

1566.000 

1566.000 

1503.521 

1616.000 

1616.000 

1504. e76 

1667.000 

1667.000 

1504.669 

1717.000 

1717. 0e0 

1505.276 

1768.000 

1768.000 

1505.919 

1616.000 

1P18.000 

1506.570 

1869.000 

1863.000 

1507.254 

1919.000 

1919.000 

1507.943 

1970  .00? 

1970.000 

1508.662 

2020 .000 

2020.000 

1509.383 

2071.000 

2071. t 00 

1510.133 

2121.000 

2121.000 

15ie.860 

2172.000 

2172.000 

1511.655 

2222.000 

2222.000 

1512.426 

2273.000 

2273.000 

1513.223 

2323.000 

2323. e00 

1514.013 
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SAMPLING 

POINTS 


1621 


HOSE  CHARACTERISTICS 


_ 

HOSE 

NUMBER 

VAU  NUMBER 

( 1/M ) 

PHASE  VELOCITT 
(M/SEC) 

CROUP  mocm 
(M/SIC) 

ATTENUATION  TERMS 

1 

.209369849 

1500.4991 

1500.0009 

.1000166031*81 

0. 

e. 

0. 

0. 

C 

.289231409 

1501.4916 

1500.0084 

. 1 8004 9429 E* 01 

0. 

0. 

8. 

0. 

3 

.209094015 

1502.4785 

1500.0227 

.1000818061*01 

0. 

0. 

0. 

0. 

4 

.208957603 

1503.4594 

1500.0443 

.1001137331*01 

0. 

0. 

e. 

0. 

5 

.208622191 

1504.4343 

1508 .0725 

.1001452231*81 

0. 

0. 

e. 

0. 

6 

.208687766 

1505.4034 

1500.1074 

.100176282E+01 

0. 

0. 

0. 

8. 

7 

.208554294 

1506.3668 

1500 .14810 

.1002069441*01 

a 

.208421764 

1507.3247 

1500.1949 

•120237195E+01 

0. 

0. 

0. 

0. 

.208290160 

1508.2770 

1500.2476 

. 1002670 54E+01 
0. 

0. 

e. 

0. 

.206169474 

1509.2240 

1500.3062 

.100296520E+01 

0. 

0. 

0. 

0. 

.208029680 

•  1510.1656 

1500.3690 

.1003256531+01 

0. 

0. 

0. 

0. 

.207900757 

1511.1021 

1500.4385 

.1003543781*01 

0. 

0. 

0. 

0. 

.207772716 

1512.0333 

1500.5118 

.  100382785E+01 
0. 

0. 

0. 

0. 

.207645512 

1512.9596 

1500.5901 

.1004108431+01 

0. 

0. 

0. 

0. 

.207519141 

1513.8609 

150e.6715 

.100438618E+01 

0. 

0. 

0. 

0. 

.207393575 

1514.9975 

1500.7568 

.1 00466094 E+ 01 

e. 

0. 


III-34 


17 

.207268802 

1515.7094 

1500.8467 

.1 00493244 E* 01 

e. 

0. 

0. 

e. 

ie 

.207144803 

1516.6167 

1500.9353 

.102520271 2*21 
0. 

0. 

0. 

0. 

19 

.207021525 

1517.5156 

1501 .0229 

.1025471902*0; 

e. 

0. 

0. 

0. 

20 

.206898936 

1518.4190 

1501 

.  1  .,05740422*01 
0  . 

0. 

0. 

e. 

21 

.206776963 

1519.3146 

1501. 1P26 

.1006010972*01 

0. 

e. 

0. 

0. 

22 

.206655527 

1520.2074 

1501.2477 

. 1£262  6377  2*01 
0. 

0. 

0. 

0. 

23 

.206534551 

1521 .0979 

1501.2999 

.1026559952*01 

0. 

0. 

0. 

0. 

24 

.206413936 

1521 .9867 

1501 .3375 

.1006040461*01 

0. 

0. 

0. 

0. 

I I 1—35 


25 

.206293595 

1522.8745 

1501.3599 

.1007125561*01 

0. 

0. 

0. 

0. 

26 

.206173442 

1523.7620 

1501.3723 

.100741378E*01 

0. 

0. 

0. 

0. 

27 

.206053429 

1524.6495 

1501 .3780 

.1007704131*01 

0. 

0. 

0. 

0. 

28 

.205933512 

1525.5374 

1501.3603 

.1007995611*01 

0. 

0. 

0. 

0. 

29 

.205813688 

1526.4255 

1501.3662 

.100326564E*01 

0. 

0. 

0. 

0. 

30 

.205693961 

1527.3140 

1501.3962 

.1008574651*01 

0. 

0. 

0. 

0. 

31 

.205574351 

1528.2026 

1501.4155 

.U0F86027E*01 

0. 

0. 

0. 

0. 

32 

.205454690 

1529.0912 

1501 .4422 

.1009143242*01 

0. 

0. 

e. 

0. 

33 

.205335596 

1529.9796 

1501 .4777 

.1009423062*01 

e. 

0. 

0. 


■  ¥ 

\ 
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34 

.205216505 

1530.8674 

1501.5261 

.1009698212*01 

0. 

0. 

0. 

e. 

35 

.205097665 

1531.7545 

1501.5867 

.1009968821*01 

0. 

e. 

0. 

0. 

36 

.204979097 

1532.6405 

1501.6534 

.  101023695E*01 
0. 

0. 

e. 

e. 

37 

.204660823 

1533.5254 

1501 .7309 

.1010500892*01 

e. 

0. 

0. 

e. 

36 

.204742872 

1534.4088 

1501.8176 

.1010761142*01 

0. 

0. 

0. 

0. 

30 

.204625265 

1535.2907 

I5ei.9117 

.101 101 821 E*01 
0. 

0. 

e. 

0. 

40 

.204508017 

1536.1709 

1502.0126 

.1011272322*01 

0. 

0. 

0. 

0. 

41 

.204391138 

1537.0403 

1502.1210 

.1011523202*01 

0. 

0. 

0. 

0. 

41 

.204274661 

1537.925 8 

1502.2427 

.1011768772*01 

Ill— 37 


r 


43 


204158598 


1 538.608 1 


1502.3648 


101201335E-01 


1 


44  .204042940 


1539.6723 


1502.4916 


45  .203927701  1540.5424  1502.6258 


46  .203812890  1541. 4ie2  1502.7650 


47  .203698514  1542.2757  1502.9104 


4e  .203584580  1543.1388  1503.0596 


49  .283471035  1543.9996  1503.2120 

i 

* 

} 

i  50  .203358035  1544.8579  1503.3699 

M 

i 

! 


\ 

i 
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0. 

e. 

e. 

e. 


e. 

0. 

0. 

0. 


.1012255601*01 

0. 

0. 

0. 

e. 


. 121 249448 E*01 
0. 
e. 

0. 

0. 


.I212?30?6E*01 

e. 

0. 

0  . 

0. 


.1012964101*21 

0. 

0. 

0. 

0. 


.1013195251*21 

0. 

0  . 

0. 

0. 


.1013424361*21 

0. 

0. 

e. 

0. 


i 


i 


.1013650721*01 
.«7H71277E-16 
0 . 

0. 

e. 


51 

.203245433 

1545.7138 

1503.5302 

.1013875351*01  ' 
.1783133581-16 
0. 

0. 

0. 

62 

.203133270 

1546.5673 

1503.6904 

.1014059061*01 

.2766399861-15 

0. 

0. 

0. 

53 

.2e3e21547 

1547.4163 

1503.8591 

.101431903E+01 

.3659426291-14 

.1416686221-19 

0. 

0. 

54 

.202910260 

1548.2669 

1504.0322 

.1014536551*01 
.436510 1971-13 
.3726873061-18 
0. 

0. 

55 

.202799455 

1549.1130 

1504.2025 

.1014754091*01 

.5597693061-12 

.623^555451-17 

C. 

0. 

56 

.202689068 

1549.9566 

1504.3776 

.Ik  149692  5t*?l 
.5669516951-11 
.7335975341-16 
0. 

0. 

5? 

.202579124 

1550.7976 

1504.5567 

.1015181771*01 
.5596171621-10 
.7566748511-15 
.46597 5273E-20 
.5283411691-20 

56 

.202469616 

1551 .6366 

1504.7334 

.1015394381+01 
.4039213461-09 
. 704550 649E-14 
.4711025061-19 
. 53442 9446E-19 

59 

.202360538 

1552.4730 

1504.9159 

.1015604081+01 

.3706642371-08 

tfimwinli  ' 


* 


111-39 


.5935011931-13 
.43121 6473E-18 
.4694330571-18 


60 

.202251899 

1553.3069 

1505.1027 

.1£1581141£*01 

.267164828E-07 

.4508160241-12 

.3560857851-17 

.4043664E6E-17 

61 

.202143694 

1554.1383 

1505.2897 

.1 016017661*01 
.169246186 E-B 6 
.3074794861-11 
.264164170E-16 
.3eei34670E-16 

62 

.202035916 

1554.6674 

1505.4796 

.1016221331*01 

.9580404562-06 

.187406181E-10 

.1752165041-15 

.19917,’122E-15 

63 

.201926572 

1555.7940 

1505.6706 

.1016420631*01 
.4815709211-05 
.U 1435324  E-09 
.1032652561-14 
.1174464811-14 

64 

.201621650 

1556.6163 

1505.8630 

.1016606031*01 
.2132295651-04 
.463651 e571-09 
.5364431561-14 
.6104223581-14 

65 

.201715140 

1557.4402 

1506.0469 

.1016748951*01 

.8224083251-34 

.202691026E-66 

.2429089851-13 

.2765482141-13 

66 

.201606995 

1558.2602 

1506.1913 

.1016775701*01 

.2719492201-03 

.7154741301-08 

.9437636231-13 

.1075005951-12 

67 

.201503035 

1559.0766 

1506.1797 

.1016562321*01 

.7543620021-03 

.2136061171-07 

.3079270761-12 

.3509261831-12 
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.201396736 

1559.9024 

1505.7967 

.1015694401+01 
.1717555961-32 
.5247444021-07 
.8262517261-12 
.9421 10374 £-12 

.261289147 

1560.7362 

1504.8423 

.101512304E+01 
•318312435E-02 
.104966052 £-66 
.  16109 5645E-11 
.2065960941-11 

.201179063 

1561.5903 

1503.2690 

.101421679E-01 
.4886953651-22 
.17431 4758E-06 
.330624690 £-11 
.3773807641-11 

.201065535 

1562.4720 

1501.4005 

. 101 356461 E+01 
.647468547 £-02 
.2505349031-06 
.5246662591-11 
.596191664E-11 

.200948140 

1563.3848 

1496.3953 

.1013268041+01 
.7751880 83E-02 
.3264496351-06 
.7585159411-11 
•566750651E-11 

.200826855 

1564.3290 

1497.4303 

.  12 132 9286 £+01 
•87e 5994991-02 
.  4 0033460 2E-26 
.103’’37566i-10 
.1x66099751-10 

.200701852 

1565.3033 

1495.5573 

.1013551601+01 
.9426034071-02 
.47381 9841 £-06 
.1376469601-12 
.157476518E-10 

.200573347 

1566.3061 

1493.7778 

.  1013S6704E+01 
.9928488001-02 
.5496065211-06 
.  1 7995410 6E- 10 
.2060075481-18 

.200441542 

1567.3361 

1492.0796 

. 12144  9USE+01 
.1033428061-01 

2339731161-10 

2680200301-10 


77 

.200306609 

1568.3919 

1460.4454 

.1015057471+01 
.1066548511-01 
.7164677401-06 
. 3 0423 4793 E- 10 
.3487365731-10 

78 

.200168687 

1569.4726 

14 ee. 86i 5 

. 101567135E-01 
.1065101351-01 
.8190871651-06 
.3970509321-10 
.4554360611-10 

79 

.200027889 

1570.5773 

1467.3167 

.1016306351+01 
.112106076E-01 
.9324267621-06 
.5213588431-10 
.59644 1023 E- 12 

60 

.199884311 

1571.7055 

1485.8921 

.1016951381+01 
.il45675ei,E-01 
.1062840831-05 
.6900044521-10 
. 7 S25c31 541-10 

81 

.199738030 

1572.8565 

1464.3112 

.1017596501+01 

.1170490321-01 

.1214231631-05 

.9217376951-10 

.1656534991-09 

82 

.199589112 

1574.0301 

1462.8381 

.1016242231+01 

.1195655061-01 

.1391222641-05 

.1244319341-09 

.143140179E-09 

83 

.196437612 

1575.2256 

14 ei. 3786 

.1018877611+01 
.1221976371-01 
.1599456611-06 
. 1699582761-09 
.1956592131-09 

84 

.109283579 

1576.4433 

1478.9301 

.1019500641+01 

.1246986591-01 

.1645908951-05 

.2351680851-06 

.2709387771-09 
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¥'■ 


85 

.199127657 

1577.6825 

1478.4895 

.1022107251*9 
.1280157121-0 
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.213357431 E-0? 
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.3607788601-01 
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.1584427551-02 

.2370810881-02 


.164623046 


.164176603 


.163724531 


.163266517 


.162602616 


.162333103 


,161858349 


161376753 


160 894714 


1908.3553 


1225.7860 


1913.5447 


1221.6244 


1918. 62e3 


1217.6648 


1924.2112 


1213.4534 


1929.6942 


1209.2873 


1936.2754 


1205.2146 


1940.9519 


1201.2527 


1946.7201 


11S7.4025 


.1173730 551+4 
.3405301941-1 
.3771269492-3 
.1566519531' 
.2356134722' 


.1 16244S57E*0 
.3110430721-0. 

.35741 iee et-@ 

.1496000561' 
.2265016851' 


.119233874E-.0J 
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.3309189931-01 
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.1202251811+0: 
.2476702451-0 
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. 1984 8304 5E-0 
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.222e?iie3r-0i| 

.2776882931-0 

.1204927751-0* 
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.1129548241-02 
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.159419964 
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217 
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220 
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.3079788422-01 
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2002.0962 
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2008.5475 
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.1231201412*01 
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.1P6U3564 
.  2  93  2  7  66441 
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.151660962 

2071.4577 

1125.2613 

.1312507162*01 

.2227601661-01 

•197619878E-01 

.1217948652-02 

.221492410F-02 

232 

.15109944? 

2079.1556 

1120 . 6960 

.1323742752*01 
.299201 647 E-01 
. 1 92 830261 2-01 
.1136764271-02 
.20-688321 1 E-02 

233 

.150532113 

2066.991? 

1116.2199 

.1327956242+01 
.200514735E-91 
.1730087S5E-01 
.10945 534 7 E-02 
.2<rl4346?9E-e2 

234 

.149959336 

2094.9630 

1111. e215 

.1334261041*01 

.iee211666E-01 

.1676043482-01 

.1056844531-02 

.1968462231-02 

235 

.149381466 

2103.0672 

1107.4944 

.1339727511*01 
.1959523491-01 
.  166277P73E-2 1 
.1056331381-02 
.2010310072-02 

236 

.I4679ee3e 

2111.3020 

1103.2378 

.  1 3443 7970 E*0i 
.1995650471-01 
.16861 6821 E-i 1 
.1083063222-02 
.2081701301-02 

23? 

.146211769 

2119.6648 

1099.0567 

. 13431 9105E*01 
.2971415611-01 
.l’75643830l-01 
. 1 13464295E-0  2 
.2207467472-02 

238 

.147620654 

2128.1525 

1094.9626 

.1 351090522*01 
.21861 5647E-01 
.1866019222-01 
.1216079932-02 
.239575966E-02 

238 

.147025903 

2136.7613 

1090.9722 

.1352968442*01 

.2349114722-01 
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2029321111-01 

133192142E-02 

2656135061-02 


240 

.146426200 

2145.4863 

1087.1257 

.1353696681.-01 
. 25573 4194E-01 
.2246959971-01 
.14S775756E-02 
.300900634  E-02 

241 

.145827492 
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.1351443401-01 
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.1936152701-02 
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.134082088 2-01 
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.2229336351-02 
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.144015772 
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.1346133711-01 

.3738734221-01 
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.2541244751-02 

.5445912241-02 
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.143429211 
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.4055627241-01 
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.6734677381-02 


247 
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2239.4712 
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248 

.141588761 

2219.1285 

1056.2721 
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.3096334041-02 
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.1375976911+01 
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250 

.140306980 

2239.0851 
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.1390437011*01 
.3006997761-01 
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.2754400791-02 
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251 

.139662538 

2249.4168 
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.1404729971+01 
.2678696031-01 
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.2236657641-00 
.5651191341-00 

254 
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.137000250 

2293.1291 

1018.4404 

.1449405681-01 
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.5667974341-00 
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.3:7'>”P560E-01 
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267 
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IV.  PLMODE:  DOCUMENTATION 


The  following  chapter  is  extracted  from  the  computer  model  PLMODE. 


PLMODE  ARL: UT  NORMAL  MODE  BASED  PROPAGATION  LOSS  MODEL 

VERSION  ).2  MARCH  i980 

oesigned  and  implemented  BY 

SUSAN  Q.  PAYNE 

DOCUMENTATION  for  plmode  IS  given  in 
the  acoustic  normal  mode  model  nemesis 

BY 

RUTH  GONZALEZ  AND  KENNETH  E.  HAWKER 
ARL: UT  TECHNICAL  REPORT  ARl-TR-80-13 
AND 

users  manual  for  nemesis  and  plmode 

BY 

RUTH  GONZALEZ  AND  SUSAw  g.  PAYNE 

arl:ut  technical  memorandum  arl-tm-bo-6 
applied  research  laboratories 
the  university  of  te*as  at  austin 
SPONSORED  BY  nORDA  code  5?0»  NSTl  station.  MISSISSIPPI  39529 


objective 


using  normal  mode  THEORY,  this  model  is  designed  to  COMPUTE  mode 
ATTENUATIONS,  complex  VELOCITY  POTENTIAL,  COHERENT  and  incoherent 

propagation  loss*  and  phase  in  a  Range "Invariant  environment. 


1 


PHYSICAL  DESCRIPTION 


PLMODE  IS  A  COMPUTER  PROGRAM  DESIGNED  FOR  USE  WITH  THE  ARL.'UT  NORMAL 
MODE  MODEL.  NEMESIS.  THIS  MODEL  ASSUMES  A  HORIZONTALLY  STRATIFIED 

deep  ocean  with  single  channel  profiles  and  multiple  fluid  sediment 
layers  overlying  a  substrate,  sound  speeds  in  the  water  and 

SEDIMENT  LAVERS  VARY  WITH  DEPTH.  THE  DEnSTTY  IS  CONSTANT  WITHIN 
each  layer,  the  last  layer  Is  a  HOMOGENEOUS,  SEMI-INFINITE 
fluid  or  solid  substrate,  i,e„  the  compress ional  and  shear 

WAVE  SPEEDS  AND  DENSITY  REMAIN  CONSTANT  WITH  DEPTH.  SOUND  SPEEDS 
ANO  DENSITIES  CAN  BE  DISCONTINUOUS  AT  LAYER  INTERFACES.  THE 

attenuations  in  this  environment  ARE  defined  IN  PLMODE,  attenuation 
IN  THE  WATER  ANO  THE  SHEAR  AND  COMPRESSlONAL  WAVE  ATTENUATIONS  IN  THE 
substrate  are  assumed  to  remain  constant  wtTh  depth,  compressional 
WAVE  ATTENUATION  £n  THE  SEoIPENT  LAYERS  CAN  yARY  WITH  DEPTH,  AND  CAN 
BE  DISCONTINUOUS  AT  THE  LAYER  INTERFACES. 
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MATHEMATICAL  DESCRIPTION 


plmode  REQUIRES  input  from  The  mode  model  tn  order  to  obtain  the 

EIGENVALUES  (KN)t  THE  EIGENFUNCTIONS  (UN(Z)),AND  THE  ATTENUATION 

terms  in  each  layer  <gnd. 

the  attenuation  terms  are  computed  by  an  integral  derived  from 
first  order  perturbation  theory  (See  arl-tp-bo-13> •  assuming  that  the 
attenuation  is  CONSTANT  in  EACH  layer,  in  PlMODE*  it  is  ASSUMED 
that  this  is  true  in  the  water  (aw,  the  attenuation  value  in  the 
WATER)  AND  The  SUBSTRATE  <AC  and  AS,  THE  COMPRESSI0NAL  and  shear  wave 
ATTENUATIONS.  RESPECTIVELY).  HOWEVER.  THE  ATTENUATIONS  MAY  VARY 
WITH  OEPTH  In  THE  sediment  LAYERS  (AL(Z)>. 

IF  THE  sediment  attenuations  are  constant  tn  each  layer  the  mode 
ATTENUATION  for  mooe  number  n  is  computed  as 

SUM  (  A  *  G  )  ♦  A  *  G  ♦  A*G  ♦  A  *  G 

L  L  Nl  w  NW  C  NC  s  nS 

where  l  Refers  to  a  sediment  layer*  w  the  water*  c  the  compressional 

WAVE  IN  THE  SUBSTRATE,  S  ThE  SHEAR  WAVE  IN  THE  SUBSTRATE,  AND  N  THE 

mode  number. 

IF  THE  SEDIMENT  ATTENUATIONS  VARY  WITH  DEPTH,  THEN  ThE  MODE 
ATTENUATION  FOR  MODE  NUMBER  Nf  ATN*  is  computed  as 

SUM  (  D  *  2  *  PI  *  F  /  K  •  INt  >  ♦  A  *  G  ♦  A  *  G  ♦  A  *  G 
L  L  N  L  W  MW  C  NC  S  NS 

2 

INT  ■  INTEGRAU  U  IZ)  •  A  (Z>  /  C(Z)  ) 

L  DZ  N  L 

WHERE  dl  IS  The  density  of  SEDIMENT  layer  l*  f  is  frequency.  c<z>  is 
the  sound  speed  as  function  of  depth,  and  tntegRal<  >dz  is  the 
integral  with  respect  to  depth. 

SINCE  ThE  ATTENUATIONS  Aw,  AC,  AS  And  AL ( z >  are  input  in  db/m/khz 
The  ATTENUATION  SUMS,  ATN*  HAVE  to  be  CONVERTED  TO  NEPERS/M 
before  use  In  THE  FIELD  CALCULATIONS* 


f 
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M  a  A  T  *  (F  /  1600.)  *  (LOG  (10.)  /  20.) 

N  N  E 

WHERE  HN  IS  the  ATTENUATION  Of  the  nth  HoDF •  AND  F  is  THE  FREQUENCY 
(IN  HERTZ). 

the  FOLLOWING  EQUATION  Is  USED  to  COMPUTE  THE  COMPLEX  VELOCITY 
POTENTIAL  (VP)  for  a  GIVEN  SOURCE  depth  <ZS),  given  RECEIVER  DEPTH 
(ZR)  AT  A  range*  R* 

I  *  PI  *  Q  *  D  *  SUH  (  H  •  u  (Z  )  *  U  (Z  )  *  EXP  (  -R  *  M  )  ) 

SNNNSNR  n 


(  SL  /  20.  ) 

Q  ■  10. 


I 

H  a  H  (  K  *  Rt 
NON 


WHERE  DS  IS  THE  OEnSITY  OF  THE  LAYER  IN  WHICH  THE  SOURCE  LIES*  SL  IS 
The  SOURCE  LEVEL  (OB  RE  1  MICROPASCAL  at  1  M).  I  IS  THE  COMPLEX  NUMBER 
l*  ftNo  H01(  )  IS  THE  HANKEL  FUNCTION  OF  ORnER  ZERO*  FIRST  KIND.  THE 

hankel  function  is  approximated  by  the  first  term  of  the  asymptotic 
expansion* 


H  (K  •  R)  *  SORT (2  /  (PI  *  K  *  R) )  *  EXP(  -I  *  (K  *  R  -  PI  /  * > >  • 

ON  N  N 

the  propagation  loss  <pl>  is  computed  by 


PL  ■  -10*  •  LOG  (  ABs(VP)  )  ♦  SL  . 

10 


The  incoherent  propagation  loss  <ipl>  is  computed  from  the  incoherent 
summation  of  the  mode  terms* 


I  *  pi  •  d  *  sum  (  (abs(h  )  *  u  (z  )  *  u  (z  )  •  exp  c  — R  ♦  m  i >  ) 

SN  NNSnR  N 


IPL  •  -10.  *  LOG  (  IMS  ) 

10 


WHERE  IMS  IS  THE  INCOHERENT  MODE  SUM. 
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The  phase  i?  computed  FROM  the  COMPLEX  Velocity  POTENTIAL* 

-1 

PHASE  «  TAN  (  IM(VP)  /  Rl(VP)  > 

WHERE  IM(VP)  AND  RL<Vp)  ARE  THE  IMAGINARY  AND  REAL  PARTS. 

respectively#  of  the  velocity  potential. 

THE  MAXIMUM  PROPAGATION  LOSS  IS  SET  AT  ISO  DB*  IF  THE  COMPLEX  FIELD 
IS  ZERO*  THE  PHASE  IS  SET  AT  10»*99. 


ON  INPUT 


there  are  two  input  piles  required  by  plmode  -  a  formatted  data 
file  supplied  by  the  USER,  and  the  disk/tape  file  createo  by  The 
normal  mode  model  NEMESIS. 


lfdisk  (Taped 


this  is  the  disk/tape  file  created  by  the  mormal  mode  model*  nemesis 
it  contains  all  the  Relevant  input  data  supplied  to  nemesis  together 
with  the  EIGENVALUES,  GROUP  VELOCITIES,  ATTENUATION  terms,  and 
NORMAL  MODES.  SEE  THE  DESCRIPTION  OF  THE  FORMAT  OF  THIS  FILE  IN 

the  nemesis  program*  or  in  subroutine  inpmod  of  plmode. 
this  file  is  read  by  routines  inpmod  and  gftmod. 


LFINP  (TAPES,  INPUT) 


This  is  The  USER-SUPPLIED  data  set  containing  the  description  of 

The  ATTENUATION  IN  the  ENVlRONMcNT,  and  THF  SOURCE/RECEIVER/RANGE 

geometries. 

this  file  is  reao  by  routines  inattn  and  inpsRr. 

the  following  restrictions  are  IMPOSED  On  the  input  to  PlMODEI 

(1)  ALL  DEPTHS  are  REFERENCED  TO  THE  AlR-WATER  INTERFACE,  AND  NONt 
CAN  fall  IN  The  SUBSTRATE. 

<2)  ATTENUATIONS  MUST  BE  DEFINED  FOR  each  Layer  that  WAS  SPECIFIED 
in  the  normal  mode  model. 

(3)  for  ATTENUATION  PROFILES,  ALL  DEPTHS  must  be  in  ORDFR  of 
INCREASING  DEPTH,  There  must  be  AT  least  A  DEPTm-aTTENUaTION 

pair  specified  on  the  upper  and  lower  interfaces  of  a  sediment 

LAYER.  THE  LAST  DEPTH  AT  ThE  BOTTOM  OF  A  LAYER  MUST  EQUAL  THE 

FIRST  DEPTH  at  The  TOP  CF  The  NEXT  LAVER. 

the  MAXIMUM  NUMBER  OF  DEPTH-ATTENUATION  PAIRS  IS  48. 

(4)  ALL  ATTENUATIONS  MUST  BE  INPUT  IN  DB/mETER/KILOHERTZ. 

<5)  ALL  DEPTHS  MUST  BE  INPUT  IN  METERS*  AMD  MUST  BE  GREATER  THAN  0* 
(6)  ALL  RANGES  MUST  BE  INPUT  IN  KILOMETERS,  AND  MUST  BE  GREATER 
THAN  0. 

<7»  THERE  ARE  NO  RESTRICTIONS  ON  THE  NUmhfR  OF  SOURCES*  RECEIVERS* 
or  Ranges  that  can  be  input  in  a  data  set  (but  see  -current 
LIMITS-  RELOW) 


input  to  plmode  is  via  formatted  read  statfments.  each  item  must  at 


POSITIONED  In  the  CORRECT  FIELD  AS  GIVEN  By  THE  FORMAT  STATEMENT. 

CARO  COLUMN  NO 

NO  OF  CARD  DESCRIPTION  OF  DATA  FORMAT 

1  1-80  header  information  used  to  describe  data  baio 

2  1-10  ATTENUATION  In  THE  water  Flo.* 

H-20  cqmpRessional  wave  attenuation  in  the  fio.* 

SUBSTRATE 

21-30  SHEAR  WAVE  ATTENUATION  In  the  SUBSTRATE  Fly.* 

IF  SEDIMENT  LAYERS  WERE  DEFINED  IN  nemesis:  CARDS  3,*  are  input* 

3  1-*  DESCRIPTION  OF  THE  TYPE  OF  ATTENUATION  InPUT*.  1A* 


CONS  TO  DENOTE  CONSTANT  ATTENUATION  IN  EACH 

SEDIMENT  LAYER* 

PROF  TO  DENOTE  DEPTh-aTtENUATIon  DESCRIPTION 

for  EACH  SEDIMENT  LAYER 

♦A.B..  FOR  CONSTANT  ATTENUATION  (CONS) * 

ATTENUATION  In  each  SEDIMENT  layer 
1-10  ATTENUATION  FOR  FIRST  SEDImEnT  LAYER 

11-20  ATTENUATION  FOR  SECOND  SEDIMENT  LAYER 

•««•  •••• 

FOR  ATTENUATION  PROFILES  (PRoF>* 
depth.aTtenuation  pairs,  there  must  be  a  pair 
defined  at  the  TOP  and  bottom  of  EACH  sediment 


LAYER 

1-10  DEPTH  Flo. 2 
11-20  SEDIMENT  ATTENUATION  fio.* 
21-30  DEPTH  Flo.2 
31-*0  SEDIMENT  ATTENUATION  Flo.* 


•••*  •••• 

5  BLANK  CARD 

an  ENd-OF-fILE  FOLLOWING  Card  5  INDICATES  the 
end  of  the  data  set  (no  somrce/Receiver/range 
input  to  COME).  CARDS  1-5-EOF  CAN  BE  REPEATED 
FOR  DIFFERENT  SETS  OF  aTTEnUaTiOnS  TO  COMPUTE 

mode  attenuations  only. 


Fly.* 

Flo.* 

Fly.* 


V-  \ 


IV-10 


7A»B. • 


1-10 

11-20 


9A.B. • 


designator  of  the  type  of  output  desired, 
depth  to  denote  FIELD,  propagation  loss 
•Vs.  depth  for  FACH  source/range 
combination 

range  TO  denote  field,  propagation  LOSS 

•VS.  RANGE  FOR  FACH  SOURCE/RECEIVER 

combination 

both  to  denote  both  types  of  output 

SOURCE  DEPTHS 
SOURCE  depth 
SOURCE  depth 


Flo. 2 
Flo. 2 


BLANK  card 

RECEIVER  DEPTHS.  8F1m.2 

can  be  input  as  individual  depths  (e.g.. 

10.  1 025, ) .  ORt 

cak  be  input  as  triplets  of  depth  increment 
(PrECEOED  by  A  NEGATIVE  SIGN).  START  DEPTH  AND 
STOP  DEPTH  (E.G.»  “10.  loo.  ?000.>*  OR» 

CAN  Be  INPUT  as  COMBINATIONS  OF  THESE  TWO  FORMS 
OF  INPUT  (E.G. .  1025.  -10.  loO.  2000.  10.  ) 

BLANK  CARD 

RANGES.  8F10.2 

can  be  Input  in  the  Forma i  described 
ABOVE  for  the  receivers. 

BLANK  CARO 

CARDS  6-12  cam  BE  REPEATED  ▼  INPUT  DIFFERENT 
SOuRCE/RECElVER/RANGE  data  >  CORRESPONDING 
TO  THE  ATTENUATION  DATA  SET  INPUT  ON  CARDS  1-5. 

AN  ENd-OF-fILE  FOLLOWING  Card  12  INDICATES  the 
eno  of  the  entire  data  set. 

THE  NEXT  CARD  TO  BE  READ  FOLLOWING  AN  END-OF-FILE 
HERE  IS  CARO  l . 


IV-11 


_ ; _ d  da 


ON  OUTPUT 


LFOUT  (TAPE6) 


this  is  the  printeo  output  file*  it  contains  the  geoacoustic  profile 
describing  the  environment,  the  mode  attenuation  table,  the  complex 
velocity  potential*  coherent  and  INCOHERENT  PROPAGATION  LOSS*  and 
PHASE  FOR  each  SOURCE/RECEIVER/RANGE  COMBINATION. 

This  file  is  written  on  by  routines  plout*  plouTi*  plOutb*  plhead* 

PLPRNT*  ALERT*  AND  ERROR. 

THE  CHARACTER  STRUGS  BOD  AND  EOD  (BEGINNING  AND  END  OF  DATA) 
DELINEATE  TH£  START  AND  ENn  OF  EACH  TABLE  "F  OUTPUT.  NOTE  THAT  WHtN 
The  phase  of  THE  ACOUSTIC  field  IS  INDETERMINATE  (COMPLEX  FIELD 
EQUALS  zero).  ITS  value  IS  SET  AT  lo##99.  ThIS  NUMBER  EXCEEDS  THE 
OUTPUT  FIELD  PROVIqEO  FOR  IN  THL  FORMAT*  AND  ON  THE  CYBER  171* 
ASTERISKS  ARE  OUTPUT. 

The  following  error  messages  may  be  printed  0n  detection  of  an  inpui 

OR  EXECUTION  PROBLEM) 


ARRAY  size  diagnostics 

insufficient  SPACE  ALLOTTED  To  ARRAY  ITITLF  in  main  PROGRAM 
the  VALUE  of  MAxTLe  In  MAIN  PROGRAM  Plmode  IS  TOO  SMALL:  IT 
MUST  be  GREATER  than  OR  EQUAL  THE  NUMflto  OF  WORDS  IN  THE 

first  record  on  the  nemesis  output  file. 

INSUFFICIENT  SPACE  ALLOTTED  To  ARRAY  BuF 

THE  DIMENSION  OF  ARRAY  rUF  IN  SUbROUTiNf  INPMOO  (NBUF)  IS 

too  small:  it  must  he  greater  Than  or  eoual  to  the  maximum 
number  of  layers  specified  in  nemesis  amd  the  number  of  words 
in  the  second  Record  on  the  nemesis  disk  file. 

MAX  NO.  UF  LAYERS  ON  MODE  DISK  .GT.  NO.  SET  FOR  THIS  PROGRAM 
THE  VALUE  OF  MAXLAY  IN  MAIN  PROGRAM  Plmode  IS  TOO  SMALL*  IT 
MUST  BE  GREATER  THaN  OR  EQUAL  TO  THE  MAXIMUM  NUMBER  OF  LAYERS 

specified  in  nemesis. 

insufficient  space  in  work  array  FOR  modes  and  attenuations 

INCREASE  The  DIMENSION  OF  kORK  ARRAY  (MaXwRK)  IN  MAIN  PROGRAM 
PLMODE.  SEE  CURRENT  LIMITS  (BELOW)  AS  A  GUIDE  TO  MAXWRK. 

MAX  NO.  OF  SVP  PTS.  On  MODE  U I S*  .GT.  NO.  SET  FOR  THIS  PROGRAM 
THE  VALUE  OF  MAXSVP  In  MAIN  PROGRAM  PlMoDE  IS  TOO  SMALL *  IT 
MUST  BE  GREATER  THAN  OR  EQUAL  TO  THE  MAXIMUM  NUMBER  OF  SOUND 
f  SPEED  PAIRS  SPECIFIED  In  NEMlSIS. 


IV-12 


INSUFFICIENT  WORK  SPACE  when  ALLOCATING  SPACE  for  atten  calcs. 

THE  VALUE  OF  MAXWRK  IN  MAIN  PROGRAM  PlMoDE  IS  TOO  SMALL S 
SEE  CURRENT  LIMITS  (BELOW)  AS  A  GUIDE  To  mAXWRK. 

Insufficient  space  in  attenuation  arrays:  must  be  .ge.  layers-2 

INCREASE  THE  DIMENSION  OF  aTTSED,  ATTdEo  ARRAYS  (MAXATT)  IN  MAIN 
PROGRAM  PLMODE  TO  AT  least  NUMBER  OF  LAVERS  -  2 

OVERFLOW  in  ATTEN  ARRAY  ON  ASSIGNING  ATTen  To  INTERFACES 

INCREASE  SIZE  OF  ATTSED  AND  ATTDEP  ARRAYS  (MAXATT)  IN  MAIN 
PROGRAM  PLMODE  TO  AT  LEAST  ?*NUM0ER  Of  SEDIMENT  LAYERS 

Insufficient  space  in  work  array  for  sedimfnt  atten  calcs 

INCREASE  the  DIMENSION  OF  THE  work  ARRAY  (MAXwRK )  In  main 
PROGRAM  PLMODE,  see  CURRENT  LIMITS  (BELOW)  AS  A  GUIDE  TO 
MAXWRK. 


INPUT  DIAGNOSTICS 

ERROR  DETECTED  WHIlE  READING  FROM  MODE  FlLF 

PROBLEM  WITH  THE  MODE  DISK  FILE:  CHECK  THaT  THE  FILE  HAS  THE 
CORRECT  STRUCTURE*  NUMBER, AND  SIZE  OF  RECORDS. 

EOF  ENCOUNTERED  WHILE  READING  WATER  AND  SUBSTRATE  ATTENUATIONS 

NEGATIVE  WATER  OR  SUBSTRATE  ATTENUATIONS  DETECTED  ON  IrjPUT 

EOF  ENCOUNTERED  WHILE  READING  SLDlMENT  ATlENUATlON  TYPE 

UNKNOWN  SEDIMENT  ATTENUATION  TYPE  DETECTED  (CONS  OH  PROF  ONLY) 

EOF  ENCOUNTERED  WHILE  READING  SEDIMENT  ATTENUATIONS 

NEGATIVE  SEDIMENT  atTENUaTIONS  DETECTED 

EOF  ENCOUNTERED  WHILE  READING  ATTENUATION  PROFILE 

zero  depth  found  in  attenuation  profile 

IF  BOTH  DEPTH  AND  ATTEN  ARF  ZERO,  THIS  MESSAGE  WILL  NOT  OCCUR  I  II 
IS  EXECUTED  ONLY  IF  DEPTH  IS  ZEHO  AND  ATTEN  IS  NON-ZERO. 

OVERFLOW  IN  ATTEN  ARRAYS;  oR»  BLANK  CARD  MUST  FOLLOW  ATTENS 

INCREASE  THE  SIZE  OF  ATTSCO,  ATTDEP  ARRAYS  (MAXATT)  IN  MAIN 
PROGRAM,  OR  REDUCE  THE  NUMBER  OF  PROFtLF  POINTS  INPUT. 

CHECK  THAT  A  BLANK  CARO  FOLLOWS  THE  ATTENUATION  INPUT 

NEGATIVE  VALUES  DETECTED  IN  ATTENUATION  PpOFILE 


